< 



A THEORY OF Li-DISSIPATIVE SOLVERS FOR SCALAR 
CONSERVATION LAWS WITH DISCONTINUOUS FLUX 

B. ANDREIANOV, K. H. KARLSEN, AND N. H. RISEBRO 



Abstract. We propose a general framework for the study of contractive 
semigroups of solutions to conservation laws with discontinuous flux: 

ut + f(x,u)^ = 0, f{x,u) = {^'^'''^' (CL) 



J'-iu), x>0, 

where the fluxes are mainly assumed to be continuous. Developing the 

ideas of a number of preceding works (Baiti and Jenssen 1141 . Audusse and 

» I Perthame 1121 . Garavello et al. 1351 , Adimurthi et al. [3], Biirger et al. I21| ). 

O I we claim that the whole admissibility issue is reduced to the selection of a 

family of "elementary solutions" , which are piecewise constant weak solutions 

^ of the form 

CO c(x) = c'll{^<0} +c''ll{^>0}- 

(N We refer to such a family as a "germ" . It is well known that (CL) admits many 

different contractive semigroups, some of which reflects different physical 
applications. We revisit a number of the existing admissibility (or entropy) 
, conditions and identify the germs that underly these conditions. We devote 

■^[^ specific attention to the "vanishing viscosity" germ, which is a way to express 

the "F-condition" of Diehl 1321 . For any given germ, we formulate "germ- 
based" admissibility conditions in the form of a trace condition on the fiux 
discontinuity line {x = 0} (in the spirit of Vol'pert I80| ) and in the form of 
G a family of global entropy inequalities (following Kruzhkov 1501 and Carrillo 

I I |22n. We characterize those germs that lead to the -contraction property 

for the associated admissible solutions. Our approach offers a streamlined and 
unifying perspective on many of the known entropy conditions, making it pos- 
^ sible to recover earlier uniqueness results under weaker conditions than before, 

and to provide new results for other less studied problems. Several strategies 
^*~^ for proving the existence of admissible solutions are discussed, and existence 

results are given for fiuxes satisfying some additional conditions. These are 
based on convergence results either for the vanishing viscosity method (with 
• standard viscosity or with specific viscosities "adapted" to the choice of a 

\l germ), or for specific germ-adapted finite volume schemes. 
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1. Introduction 

1.1. L^-dissipativity for scalar conservation laws. We are interested in the 
well-posedness of the Cauchy problem for a general scalar conservation law: 

ut + div f{t, X, u) =0, (1.1) 

u\t=a = uq. (1.2) 

It is well known that classical solutions to this problem may not exist globally since 



discontinuities can develop in finite time; hence (1.1 1,(1.2) must be interpreted 
in the weak (distributional) sense. However, weak solutions are in general not 
unique, and so additional admissibility criteria are needed to single out a unique 
solution; these so-called entropy conditions are usually motivated by a very careful 
inspection of the underlying physical phenomena in presence of discontinuous or 
rapidly changing solutions. 

Consider solutions which take values in a closed interval U C R. Whenever the 
flux function f : M+ x x J7 i— >■ is sufhcicntly regular in all variables, the clas- 
sical notion of Kruzhkov entropy solutions in the L°° framework [SD] (cf. Volpert 
[80] for the BV setting and Olemik fOT ) and its further extensions to bounded do- 
mains, measure-valued solutions, renormalized solutions, and degenerate parabolic 
problems, cf., e.g., HSl HZl HH IMl ISHl ISH] , have provided a rather complete theory 
of well-posedness. An equivalent notion of kinetic solutions was later formulated in 
[56L I71j , which allowed for a deeper study of regularity and compactness properties 
of admissible solutions. By far the most studied case is the autonomous equation 

wt -I- div f (u) = (1.3) 

Under mild regularity assumptions on the flux f, it is well known that the Kruzhkov 



entropy solutions of (1.3), (1.2) (with mq € n L°°) form an L^-contractive and 
order-preserving semigroup onL^nL~[50l[72l[16l[27l[5ll[52l[l8l[3[58]. Moreover, 
this is the unique semigroup of this kind which admits all the trivial constant 



solutions of (1.3). This semigroup is generally viewed as the one representing the 
physically relevant solutions to scalar conservation laws. Solution semigroups that 
are Lipschitz in but not contractive may exist and are of physical interest (see 
|54j ) ■ but their study will not be addressed in the present paper. 

A comparable theory for conservation laws with discontinuous flux is still not 
available, although these equations have received intense attention in last fifteen 
years; see [3 1 13 [II [11 [P, ^T^, 1^ IIS [11 113 [Ml S El 131 133 133 
|151llill71llHllia[SlEni[B2, 69,7l[731[7Sl[77] (and additional references therein) 
for a number of different admissibility criteria, existence and/or uniqueness results, 
which we partially revisit in Section [4| Recently, it was pointed out explicitly by 
Adimurthi, Mishra, and Veerappa Gowda in |3] that for the case f = f{x, u) with f 
piecewise constant in x, there may exist many different L^-contractive semigroups of 
solutions to (jOJ. Different semigroups correspond to different physical phenomena 
modeled by the same equation, but with different dissipative processes occurring 
on the discontinuities of the solutions. We refer to [20] for a comparison of a 
clarifier-thickener model and a porous medium model, which leads to the same 
formal conservation law, but with two distinct semigroups of physically relevant 
solutions. Notice that even for the classical equation ( [1.3[ ), non-Kruzhkov L^- 
contractive semigroups can also be constructed (see [26l I2), by an approach much 
similar to the one of [1 [21] and the present paper. 

In this paper we formulate a streamlined, unifying framework encompassing the 



different notions of entropy solutions for (1.1). For the multi-dimensional problem 



examined in the sequel paper [3] we consider Caratheodory functions f{t, x, u) that 
are piecewise Lipschitz continuous in (t, x) and locally Lipschitz continuous in u. 
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The present paper focuses on the model one-dimensional case with 

f = fix, u) = /(u)U{,<o} + r (u)ll{,>o}, (1.4) 

and : U — > M being merely continuous (in the existing literature the fluxes 

are assumed to Lipschitz continuous). We note that a generalization to a piecewise 
constant in x and continuous in u flux function f(x, u) is straightforward. 

The initial datum uq is assumed to belong to i°°(IR; U). The presence of a source 



term s{t,x) in the conservation law (1.1) is easy to take into account, for example 
in the case C/ = M and s{t, •) e L°°(R^) for a.e. t > with \\s{t, OHl- dt < oo 
for all T > 0. For the sake of simplicity, we always take s = 0. 

Definition 1.1. A semigroup {St)t>o, St '■ D L°^{R;U), defined on a subset 



D C i°°(M; U), is called an L^-dissipative solver for (1.11,(1.4) if 
• for all uq g D, the trajectory u{t,x) := St{uQ){x) of the semigroup St gives a 
solutioT^to the problem ( |1.1[ ) with flux (1.4) and initial data Uq. 



• for all Uq, Wo £ D, u := 5t(uo) and u := St{uo) satisfy the "Kato inequality" 
Ve e 2?([0,oo) X M), ^ > 0, 

- / / {\u-u\S,t+s\g\\{u-u){f{x,u)-f{x,u))i,x]< / |uo - uol ^(0,a;). 

JR+ Jm. Js. 

(1.5) 

Letting ^ go to the characteristic function of {0,t) x M (cf., e.g., [111 IS]), we 
deduce from (1.5) the L^-contractivity property 

whenever I uo — Uq I G .^^(K), / \u — iL\{t)< / |uo — ^iol for a.e. t > 0. (1.6) 

Js. Jr 

We could also have required an L^-dissipative solver to be order-preserving, in 



which case sign(M — u) and \u — u| in ( 1.5 ),( 1.6 ) are replaced by sign^(u — u) and 
(u — u)+, respectively. It turns out that the L^-dissipative solvers that we consider 
are automatically order-preserving (cf. [28J). 

Notice that in several spatial dimensions, mere continuity of the flux f in u is 
not sufficient to obtain contractivity from the Kato inequality; see in particular 
the counterexample of Kruzhkov and Panov (55] ; cf. [111 EH [THl H [SS] for some 



sufficient conditions for (1.6) to hold 



1.2. Analysis in terms of Riemann solvers. Our primary goal is to formulate 
a convenient unifying framework for the study of different i^-dissipative solvers. 
We will work with the basic one-dimensional model equation 

ut+f{x,u)x = 0, f = f(2^, w) = /'('")ll{i.<0} + /''(w)ll{2:>0}, (1.7) 

which the existing literature has targeted as the main example for understanding 
the admissibility issue in the case of a spatially discontinuous flux. 

With the admissibility of solutions in the regions {x > 0} and {x < 0} being 
understood in the sense of Kruzhkov's entropy solutions, it remains to define the 
admissibility of a solution at the discontinuity set S = (0, oo) x {0}. As suggested by 
Garavello, Natalini, Piccoli, and Terracina [35], the admissibility issue should reduce 
to the choice of a Riemann solver at x = 0. That is, to each pair (m_, m+) G U x U 



we have to assign a weak solution u TZS{u^, u+) of (1.7) with the initial datum 



Ufi{x) ^ \ (1.8) 
\u+, X > 0. 



^i.e., a weak solution which is admissible in a sense specified later. 
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If (1.5) and (1.6) turn out to be true, such a solution should be unique. If we 
assume that 



the notion of solution for (1.7) is invariant 
under the scaling (t,x) i— {kt,kx), fc > 0, 



(Al) 



it follows that the unique solution to (1.7), (1.8) is self-similar, i.e., it only depends 
on the ratio ^ := x/t. Furthermore, (1.5) should hold with u,u being solutions 
72.iS(u_, TZS(u-,u^) of two different Riemann problems. 
If, in addition, we assume 



any solution of (1.7) is a Kruzhkov entropy solution in x (M \ {0}), (A2) 

it follows that the solution u = TZS{u-,u+) is monotone in ^ on each side of the 
discontinuity; thus, it has strong pointwise left and right traces on E = (0, oo) x {0}. 
Let us denote these traces by 7'u and j'^u, respectively. If m'''' :— 7''''^, the weak 
formulation of (1.7 1 yields the Rankine-Hugoniot condition 

fin') 



(1.9) 



Hence, we can associate to the pair (u'',u^) the following stationary weak solution 
of (|L7): 

(„,l „ ^ n 

(1.10) 




Any solution of the form ( 1.10 ), with (m', u^) G U xU, will be called an elementary 
solution" , and be identified with the pair (u^u^). 

Under assumptions (Al) and (A2), we conclude that the solution 7?.5(u_,u+) 
of the Riemann problem should consist of the following ingredients: 

• the standard Kruzhkov self-similar solution joining the state u_ 
at t = 0, X < to some state u' at t > 0, a: = 0~ ; 
the jump joining the state at a; = 0^ to a state u'' at a; = 0^ 



(RPb-sol.) 



such that (1.9) holds; 

the standard Kruzhkov self-similar solution joining the state 
at i > 0, a; = 0+ to the state u+ at t = 0, x > 0. 



Remark 1.2. In (RPb-sol.), it is convenient to consider the situation in which the 
wave fans joining m_ to u' , resp. to might contain zero-speed shocks. In this 
case, the same a.e. defined solution may correspond to two or more different pairs 
{u\u^) of intermediate states at x — 0*, and 7'''"^ may differ from u''*". 

Whenever solutions of Riemann problems are seen as trajectories of an L^- 
dissipative solver for ( |1.7[ ), the different "elementary solution" pairs {u'',u^), {u'',u^) 
that can be used in (RPb-sol) ought to satisfy 

q'iu',u'):=sign{u'~u'){f'{u')~f{u')) 

> sign(u'- - umriun - riul) =■■ q''{u^,ul. 



(1.11) 



Indeed, denote by u'''' the left and right traces on E of m = TZS{u-, u+) and let the 
test function ^ in the Kato inequality (1.5 1 converge to the characteristic function 
of (0, T) X {0}. We then easily arrive at ( |1.11 ). 

1.3. Admissibility germs and uniqueness. The property in ( |1.11[ ) of allowed 
jumps across E has been recognized in many works as crucial for the uniqueness 
of admissible solutions; it was the key ingredient of the uniqueness proofs in for 
example the works [31 1351 [H] . 

Our angle of attack is to "axiomatize" property ( |1.11[ ) . More precisely, we claim 
that admissibility can be defined directly from the choice of a set Q C U x U CM."^ 
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such that ([L9| and hold for all e g. Such a set g will be 



called an -dissipative admissibility germ (an L^D germ for short). If an L^D 



germ g admits a unique maximal extension, still satisfying (1.9) and (1.11), the 
germ is said to be definite. 

Our results in Section [3] can be formulated as follows: if is a definite germ, it 
corresponds to a unique L^-dissipative solver Sf for (1.7) such that 



(Sf-sol) 



for any trajectory u of Sf , the restrictions of u to the domains 
R+ X (—00,0) and x (0,oo) are entropy solutions 
in the sense of Kruzhkov; 

roughly speaking, the pair of left and right traces 
(7'u,7'"m) of u on x {0} belongs to g* pointwise, 
where g* is the unique maximal L^D extension of g. 



The set g is exactly the set of "elementary solutions" contained within Sf . 
Without non-degeneracy assumptions on the fluxes /' and /'', the second statement 
in (5'f-sol.) should be made precise. Following the idea of [ZlllIS], the pair of (weak) 
traces (7'u, j'^u) can be defined as a Young measure on R^, which is the "nonlinear 
weak-^ limit" (see [331131]) of [u{t,—hn),u{t,hn)), with (/i„)„>i, /i„ i 0, being a 
sequence of Lebesgue points of the map x (u(-, —x),u{-,x)). What we mean is 
that the support of this Young measure is contained in g, for a.e. t > 0. A more 
practical way to understand the boundary trace issue is suggested in Definition 
[3^ii) in Section [s] 

Proposition 3.10[ i) in Section [s] shows that the elementary solutions on which 



the admissibility criterion is based actually belong to Sf , i.e., they are admissibl^ 
From this viewpoint, our approach generalizes the original one of Kruzhkov by 
defining admissibility in terms of the local L^-dissipativity property (|L5|) with 
respect to a "small" set of solutions that are judged admissible a pnor j [ In the 
Kruzhkov case, these are the constant solutions of Ut+fiu)^ = (cf. Section|42|). In 



our context, these are the elementary solutions of the form ( 1.10 ) with {u\ u^) G t/, 
and the constant solutions in the domains M+ x (— oo,0) and M+ x (0,oo). 
Our perspective has the following practical implications: 

First, there exist in the literature a number of different admissibility criteria 
devised for particular forms of the fluxes /',/''. For some of them, the uniqueness 
and/or existence issues remained open; and it is often difficult to judge whether the 
same criteria can be useful for fluxes more general than those initially considered. 
On the other hand, if we manage to find some of the elementary solutions allowed 
by a given admissibility criterion in a given configuration of fiuxes /'^, it becomes 
possible to determine if the criterion is adequate. We give examples in Section |4j 

Second, some important admissibility criteria derive from certain regularized 
problems. For example, for the classical vanishing viscosity method, which consists 



in adding eAu to the right-hand side of (1.1 1 (cf. Hopf [H] and Kruzhkov [20]); a 



variant of the Kato inequality (1.5 1 is well known for each e > 0. Hence, the Kato 
inequality is easily passed on to the limiting problem, which in turn gives raise to 
an implicitly defined L^-dissipative solver for the limit problem. Determining the 



■^Notice that if Q is not L^D, then the admissibility of some of its elementary solutions is 
contradictory: the pairs (u',u^), {u'-,u^) for which fails rule out each other! 

''a well-known alternative (see Benilan 1161 and Crandall 27 ) is to define admissibility by 
the global L^-contraction property with respect to a dense set of solutions of the stationary 

problem u -f f{u)x = s; this is the nonlinear semigroup approach to scalar conservation laws. We 
will not develop this second possibility here. From another point of view, however, the definition 
of admissibility on the flux discontinuity line {x = 0} by property is a direct generalization 

of the admissibility approach of Vol'pert |80| : see Section 4.2 
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elementary solutions (as £ — > limits of standing-wave profiles, for example) opens 
up for the possibility of characterizing the germ of the associated L'^-dissipative 
solver, and as a consequence describe it explicitly. This can eventually lead to an 
existence result. Examples are given in Sections [5] and |6] Our analysis utilizes 
several simple properties of germs listed in Sections |3] and |4] 

Finally, we point out that while the formulation (S'f-sol.) (cf. Definition 3 



convenient for the uniqueness proof, it is rather ill-suited for passaging to the limit 
in a sequence of approximate solutions. Therefore we supply another formulation 
in terms of global entropy inequalities, inspired by the previous works of Baiti and 
Jenssen [14J, Audusse and Perthame [I2j, Biirger, Karlsen, and Towers and by 
the founding papers of Kruzhkov [SD] , Otto [53] , and Carrillo [52] . This formulation 



(cf. Definition 3.15 for the precise notion) is also based upon a predefined L^D germ 
Q associated with the flux discontinuity. Stability of this formulation with respect 
to Ll^^ convergence of sequences of (approximate) solutions makes it convenient for 
the existence analysis. It is shown in Theorem |3.18| that the two formulations are 
equivalent. We call the associated solutions ^/-entropy solutions. 



1.4. Measure- valued t?-entropy solutions and convergence results. The 

definition of ^/-entropy solutions in terms of global entropy inequalities can be 
adapted to define measure-valued ^/-entropy solutions (more precisely, we use the 
equivalent device of process solutions developed in [33', '34^ ; the same idea appeared 
previously in |65| ) . Although we are not able to conceive a direct uniqueness proof 
for fj-entropy process solutions, in Section[3]we prove their uniqueness if we already 
know the existence of C/-entropy solutions. A result like this makes it possible to 
prove strong convergence of merely bounded sequences of approximate ^/-entropy 
solutions. Heuristically, on the condition that existence can be established by some 
approximation method enjoying strong compactness properties (cf. Section 6.3), 
we can deduce strong convergence of other approximation methods with mere L°° 
estimates. We give the details in Section [6^ 



1.5. Outline of the paper. In Section [2| we recall some parts of the Kruzhkov 
theory; we also state and slightly reformulate some results by Panov on initial and 
boundary traces of Kruzhkov entropy solutions. In Section [3| we define admissi- 
bility germs and related notions, provide several definitions of ^/-entropy solutions, 
discuss the relations between them, provide uniqueness, -contraction, and com- 
parison results, and L°° estimates. In Section [4| we derive a series of important 
properties of germs. Moreover, we classify several known admissibility criteria in 
terms of their underlying germs, and determine their areas of applicability. An 
important example is the "vanishing viscosity" admissibility condition. In Section 
[5] we characterize its germ for general flux functions /'■''. In Section [sj we justify 
the convergence of the standard vanishing viscosity approximations to the "van- 
ishing viscosity germ" entropy solutions obtained in Section [5| Furthermore, we 
discuss variants of the vanishing viscosity approach that lead to different germs; 
in particular, we give a simplified proof of existence for the "(A, i?)-connection" 
admissibility criteria studied by Adimurthi, Mishra, Veerappa Gowda in [3] and by 
Biirger, Karlsen and Towers in [21]. Then we give an existence result for ^-entropy 
solutions corresponding to a general germ. The proof is a by-product of a strong 
convergence result for a particular germ-preserving finite volume scheme. Finally, 
we discuss convergence of L°° bounded sequences of approximate solutions towards 
C/-entropy process solutions. 
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2. Preliminaries 

First, we recall a localized version of the definition due to Kruzhkov [50]. Let 
17 C M+ X M be open and set d°n := dnn{t = 0}. For / e C{U; M; ), a function 
u E L°°(ri; U) satisfying the inequalities 

VfceC/ Vee2?(M+xM), e>0, ?bn\aoa = 0, 

// {|w-fc|6+g(u,fc)e:.} - // sign(M-fc)sC> / |uo-fc|e, 
JJn Jjr2 Jd°n 

is called Kruzhkov entropy solution of the conservation law 

ut + f{u)x^s in 17, u\ooQ = uq. (2.1) 
Here for fc G M, the function z i— > |z — fc| is called a Kruzhkov entropy, and 

q{z, k) := sign(u - k){f{u) - /(fc)) 

is the associated entropy flux. 

An important ingredient of our formulation is the fact that entropy solutions 
admit strong traces, in an appropriate sense. For example, when il = x M, an 



entropy solution of (2.1 ) admits uq as the initial trace in the sense 



VCeX'(M) lim^^// £,[x)\u{t,x) -ua{x)\^Q. (2.2) 

'(0,/i)xI 



/i^o+ h 



(see Panov [ST] and the previous works of Chen, Rascle [25 and of Vasseur [75]). 
Thus t = can be seen as a Lebesgue point of the map t ^ u{t, ■) e Ll^^{n). 

Most importantly for the present paper, analogous strong trace results hold for 
entropy solutions in half-space domains. For example, with 12 = x (— oo,0), 
provided the non-degeneracy assumption 

for any non-empty interval (a, fe) C U, f\(a,b) is not constant (2.3) 

holds, there exists a measurable function t H> (7'u)(i) on M+ such that 

Vee2?(0,oo) lim yff £_{t) \u{t,x)- {j^u){t)\^0. (2.4) 

This result is a particular case of Panov |68| Theorem 1.4] (see also [79j l53]). 
Observe that we require from the flux / only that it is continuous and satisfies the 



non-degeneracy condition ( |2.3[ ). The function 7'u is the strong left trace of the 
entropy solution u on {x — 0}. 

We emphasize that (|2.4| is a sufficient — but not a necessary — property for 



our purposes. Actually, strong traces of f{u) and of the entropy fluxes g(u, k) exist 



even in the absence of the non-degeneracy assumption (2.3). To state the result, 
we consider a "singular mapping" : [/ i-> M. There are different choices of the 
singular mappings; in particular, whenever / is of bounded variation it is convenient 
to define ^'(z) — \df{s)\ (cf. Temple [75], Klingenberg and Risebro and also 
[771 m [3]). Let us fix the choice 

Viz) := f llB(s)ds, (2.5) 
Jo 

where E is the maximal subset of U such that /(•) is non-constant on any non- 
degenerate interval contained in E. Clearly, under the non-degeneracy assumption 



(2.3), V becomes the identity mapping. Now the result of [551 Theorem 1.4] can be 



restated in the following way. 
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Theorem 2.1. Let f : U ^ R be continuous. Assume u is a Kruzhkov entropy 
solution of Ut + f{u)x = in := M-|_ x (— oo,0). Then V{u) admits a strong left 
trace on {x — 0}. Namely, there exists a function t H- {'y''V(u))(t) such that 



veel?(0,oo), 



lim — 

h^0+ h 



Similar statements hold with 



at)\V{u){t,x)-{j'V{um\^0. (2.6) 

+ x(-h,0) 

= X (0,oo) and the right trace operator 7'". 



Now let us turn to our problem ( |1.7[ ) . Whenever u is a Kruzhkov entropy solution 
of (1.71 in the domains {x < 0} and {x > 0}, by Theorem 2.1 there exist strong 
one-sided traces '~f^V^{u), 7''y(u). Here the singular mappings V''^ are defined by 
(2.5) from / = /'''', respectively. 

In the two subsequent remarks, we indicate how to use the strong traces of V'''^{u) 
to express the strong traces of the fluxes f'''^{u) and of the Kruzhkov entropy fluxes 
q^-''-{u,k) (cf. TTl). 

Remark 2.2. Define the monotone multivalued functions \V^]~^ ,\V^]~^ . Notice 
that the superpositions g^'''~ :— o are continuous functions. Moreover, 

the entropy flux functions q'"'^ , naturally written in terms of the unknown z and of 
a parameter k, can be expressed as continuous functions of V'' '"^ (z) , '"^ (k) only. 
More precisely, we have 

q\z,k)=sign{z-k) {f\z)-f\k)) 

= sign(y'(z) - V\k)) [g' o V\z) - g' o V\k)) 
=: Q\V\z),V\k)); 

a corresponding representation ofq'^{-,-) by Q^{V^{-),V^{-)) is valid. 



Remark 2.3. It is easily seen that the traces in the sense of (2.6) can be composed 
by any continuous function; namely, j'" '^ {h {V{u)) = h{'j^''^V{u)) for h G C{U; 



It follows that in the context of Theorem 2.1 , there exist strong traces ^^'^q^''^{u, k) of 



the entropy fluxes q '^{u,k), and these traces are equal to Q '^{j '^V ''^{u),V '^{k)) 
This remains true if we replace k by another entropy solution u. In this case, 

-f''''q'''\u,u) = Q'''-(V''-l^'''-(u),7'''-l/'^'-(u)). 



Recall that whenever /'''" are non- degenerate in the sense of (2.3 1, V'"'^ are just 
the identity mappings and Q'"^ coincide with g''^. Finally, notice that the Rankine- 
Hugoniot relation for a weak solution u of (1.7) can be expressed under the form 
g\-i'V\u))=g-{YV^{u)). 



3. The model one-dimensional problem 

Consider problem ( |1.7| , \1.2\ . Here the flux f is constant in x on each side of 
the discontinuity line E := (0,cx)) x {0}, on which the admissibility condition will 
be defined. As in Section [2j we denote the Kruzhkov entropy flux by 

q(x, z, k) :— sign(z — k){^{x, z) — f(a:, k)). (3.1) 

We write q^^^{z,k) for the left and right entropy fluxes sign(z — fc)(/'''"(2:) — /''''(fc)). 
Whenever it is convenient, we represent g'''' by means of the continuous func- 



2.2 



tions g^'^ , Q'''' and the singular mappings V''''^ introduced in Remark 

/'^'■(z) = .g'''^(y'''-z), g'^''(z, k) = Q^'^^iV^'^-z, V^'^'k), (3.2) 
where to simplify the notation we write V''z instead of V^{z), et cetera. 
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3.1. Definitions of germs and their basic properties. Related to the left and 
right fluxes /' and /'', we introduce the following definitions. 

Definition 3.1. Any set Q of pairs {v} ^u^) S UxU satisfying the Rankine-Hugoniot 



relation (1.91 is called an admissibility germ (a germ for short). If, in addition, 
( |1.11[ ) holds for all {u\ u^), {u\ ii^) G Q , then the germ Q is called an L^-dissipative 
admissibility germ (an L^D germ for short). 

In the sequel, we focus on L^D germs, which are those leading to the L^- 
dissipativity properties ( |1.5[ ), ( |1.6[ ). In this case, each pair {u'',u^) G Q corresponds 



to a solution of (1.7) of the form (1.10) which will be judged admissible a priori 
(see also Proposition 3.10[ i) below). 



According to the analysis carried out in the introduction, if (u', u^) G Q is judged 
to be an admissible elementary solution on (1.7), then inequality ( |1.11[ ) together 



with the Rankine-Hugoniot relation should hold for any other admissible elementary 
solution {u\u^) of (1.7). Therefore we introduce the following definition: 



Definition 3.2. Let Q be a germ. The dual germ of Q, denoted by Q* , is the 
set of pairs {u\u^) £ U x U such that ( |1.11[ ) holds for all {u\u^) G Q, and the 
Rankine-Hugoniot relation f\'u!') = f^{u^) is satisfied. 



Hence, each elementary solution of (1.7), whic h is e xpected to be admissible, 
corresponds to a pair g Q* (cf. Proposition 3.10 'ii) below). 

If Qi^ Q2 are two germs such that Qi C Q2, we say that Q2 is an extension of Qi. 
If both Qi, Q2 are L^D germs, we call Q2 an L^D extension of Qi. 

Definition 3.3. // Q is an L^D germ which does not possess a nontrivial L^D 
extension, then Q is called a maximal L^D germ. If Q is a germ that possesses a 
unique maximal L^D extension, then Q is called a definite germ. 

Notice that any maximal L^D germ is definite. As the following proposition 
shows, the definiteness of Q is necessary and sufficient for its dual germ Q* to be 
an L^D germ. 

Proposition 3.4. Fix a germ Q, and let Q* be the dual germ of Q . 

(i) One has Q <Z Q* if and only if Q is an L^D germ. 

(ii) Assume Q is an L^D germ. Then Q* is the union of all L^D extensions of Q . 
Specifically, Q is a definite germ implies Q* is Q 's unique maximal L^D extension. 
(Hi) One has Q* = Q if and only if G is a maximal L^D germ. 

(iv) If Q is a definite germ, then [Q*)* — Q* . 

(v) If Q* is an L^D germ, then Q is definite. 

Proof. Property (i) follows directly from the definitions. For a proof of (ii), let 
Q' be an L^D extension of Q. Clearly, Q' C G* ■ Reciprocally, let (u^iT) e Q* . 
Then Q' := Q \J {(«',«'")} is an L^D extension of Q which contains {u\vJ'). The 
second part of the assertion (ii) follows immediately. Property (iii) follows from (ii). 
Properties (ii) and (iii) imply (iv). For a proof of (v), we reason by contradiction. 
Let Qi, Q2 be two different maximal L^D extensions of Q. Then Q1VJQ2 is not L^D. 
By (ii), Q* contains Qi U Q2 and therefore it is not an L^D germ either. □ 

Let , y be the singular mappings introduced in Section [2) By Ran(y'''') we 



mean the images of U by the functions F'''', respectively. It is clear from ( |3.2[ ) that 
the validity of the inequality ( |1.11[ ) and of the Rankine-Hugoniot condition only 
depends on u\v},u^,ii^ through the corresponding values V^v} ^V^v} ^V^u^ ,V'''vJ' 
This motivates the following definition. 

Definition 3.5. Let Q be a germ. The reduced germ of Q is the set 

vg -.^ {{v'u^v'-u'-) \ {u\u'-)Gg]. 
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Remark 3.6. The following properties are easily derived: 

(i) If Q is a maximal L^D germ, then Q ~ | (F'u', V^'-it' ) £ VQ^. 

(a) Ifg* is the dual of a germ Q, then Q* = {{u^u') \ {V^u^V^u") e VQ*). 



Remark 3.7. In accordance with (3.2), VQ* can he equivalently defined as the set 
of pairs {{>'', v^) S Raii(y') x Ran(V^'') such that the Rankine-Hugoniot condition 
g'-{v^) —g^[v^) holds, and 

Q\v\v^) > Q''K,r), whenever {v\v'') G VQ. (3.3) 

Similarly, if Q* is anL^D germ, then for all (w', w'' ),(«', u'') S VG* inequality (3.3) 



is satisfied. Furthermore, if Q* is a maximal L^D germ, then VG* does not possess 



a nontrivial extension V' C Ran(V^ ) x Raii(y") such that (3.3) still holds for all 



In the sequel, the reader may assume that are not constant on any nontrivial 
interval, so that in all subsequent statements we have V''"'^ = Id, Q'"'" = g'''', and 
the reduced germs VG, VG* can be replaced with G, G* , respectively. 

3.2. Definitions and uniqueness of Cy-entropy solutions. We are now in a 



position to define fj-entropy solutions of (1.7), (1.2) and study their uniqueness. 

Definition 3.8. Let G be an L^D germ, with dual germ G* ■ A function u{t,x) in 
L°°(M+ X M; U) is called a G-entropy solution of (fO) , (fOl) if: 



(i) the restriction of u to il' := M+ x (—00,0) is a Kruzhkov entropy solution of 
the conservation law with flux /'(■); the restriction of u to 57'' := M_|_ x (0, oo) is a 
Kruzhkov entropy solution of the conservation law with flux /' (•); 

(ii) H^-a.e. onY. = (0, oo) x {0}, the pair of strong traces (7'V^'(m), 7''y (u)) on 
S belongs to the reduced germ VG* of G* ; 

(Hi) TL^-a.e. on {0} x M, the trace ^^u equals uq. 

The above definition fits the expectations of the preliminary analysis carried out 
in Section[l.2| Indeed, notice the following points that we list in a remark. 



Remark 3.9. 

(i) The existence of strong traces "f^V^u, j^V^u, j^u in Definition 3.8{ii) is not 
a restriction: due to Theorem 2.1 and (2.2), it follows from the point (i) of the 



same definition. Moreover, Definition 3.8 (i) and the result of i67j imply that u has 
a representative in C([0, oo); ^[^^(Ili)). In the sequel, we mean that a solution u is 
defined for all t as an L°° (]R) function. 

(ii) A G-cntropy solution of ( |1.7[ ),( 1.2 ) is a weak (distributional) solution. 

(Hi) A G-entropy solution possesses both the scaling invariance property (Al) and 



property (A2) required in Section 1.2 



(iv) In view of Remark \3.6^ ii) , any function of the form (RPb-sol.) is a G-entropy 
solution of the Riemann problem (1.7), (1.8) if and only if {u^u^) € G* . 

(v) One should compare Definition \3.8Y ii) with the definition of admissibility of 
jumps of BV solutions introduced by Vol 'pert in [80] , in the special case /' = 



(see inequality (4.1 1 in Section 4.2). 



For a proof of Remark 3.9 



ii), we notice that Kruzhkov entropy solutions are 
weak solutions on their domains; thus u is a weak solution on (0, oo) x (IR\ {0}). In 
addition, the Rankine-Hugoniot relation o n £ = M_|_ x {0} holds, because we have 

YV'u) e VG*, and Definitions 
[2?^ [2^ imply that 



3.2 



3.5 



of VG* together with Remarks 
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for some {w\w^) € Q*; in particular, f{w'-) — P'{w^). 

From Remark iv) and Proposition |3 .4| [i) , we directly deduce that the germ Q 
(and, more generally, Q*) indeed corresponds to a selection of admissible elementary 



solutions (1.10 1, as we point out in the succeeding proposition. 



then the junction u :— u^^{x<a} + u'^ixX)} Q-entropy 



of the above form is a Q-entropy solution of (1.7 1 if and only 



Proposition 3.10. 

(t) lf{u\un e g, 

solution of ( |1.7[ ). 
(ii) A function u{x) 

if {u\u'') e g*. 

Consider the map which associates to uo G L°°{U) the value u{t, •) on the 
trajectory of some C/-entropy solution u defined globally in time. In case such a 
solution exists, we state that uq belongs to D, the domain of Sf . Provided the 
underlying germ g is definite, the following theorem shows that such a solution u 
is necessarily unique; one then deduces that the domain of Sf is independent of t, 
that (Sf ) t>o is a semigroup, and it is strongly continuous in Lj^^(M), according to 
Definition 3.8 iii). In fact, the theorem shows that the map Sf is an i^-dissipative 
solver in the sense of Definition 11.11 



Theorem 3.11. Assume that g is a definite germ. If u and u are two g-entropy 



solutions of problem (1.7 1,(1.2 1 corresponding to the initial data uq and uq respec- 



tively, then the Kato inequality (1.5) and the L^-contractivity property (1.6) hold. 



In particular, there exists at most one g-entropy solution of problem (1.7), (1.2) 



Proof. We only prove the Kato inequality; the L^-contractivity ( 1.6 ) and the unique- 
ness will follow by considering in (1.5) the test functions £,fi := min{l, (i? — |a;|)"'"} 
with R — > 00, using the continuity of /'' and the fact that the space of x has 
dimension one (see [Ml [51]). 

Take ^ S 2?([0, 00) x K), ^ > 0. By a standard approximation argument, for 



} in the Kruzhkov 



h > we can take the test function — ^ min{l, ^ 

entropy formulation, in each of the subdomains il'''". By the "doubling-of- variables" 
argument of Kruzhkov [50] , we obtain the standard Kato inequality 

{\u~u\{£_h)t + q{x,u,u}{S,h)x} - / \uQ-uo\^h{0,x) <0. 



Clearly, £^h, {£,h)t converge to ^, ^4, respectively, in iy-'^(M_|- x M). Using Remark 
and calculating {^h)x explicitly, with the Landau notation Oh^o, we deduce 



2.2 



1 
h 



Jr 

Q\v'{u)y{um 

+ x(-2h,-h) 

Iff Q^-{r-{u),v^{u))^<o. 



(3.4) 



S+x(h,2h) 

Sending /i — >■ in the latter two terms, keeping in mind the definition of strong 
traces of V^''^{u),V''''^{u) and bringing into service the continuity of Q'''' (cf. (2.6) 
and Remark 2.3), we obtain 



{Q\v\t),v\t))-Q^{v^t),v%t)))at,0), 



(3.5) 



where v'' — -0' — 7'^'?)', 

sense for a.e. t > 0. By Definition 
Now, because g is assumed to be definite, by Proposition |3.4| it follows that g* is 



3.8 



— ^'"V^u^, and — ^^V^u^ in the pointwise 
"ii), we have {v\t),v'- (t)),(v\t),v'^ {t)) G Vg* . 
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an L^D germ. By Remark 3.7 we conclude that the term in (3.51 is nonnegative 
Therefore the Kato inequahty (1.51 foUows from (3.4 1 at the hmit /i — > 0. 



□ 



We collect a few comments in the next two remarks. 



Remark 3.12. 

(i) The assumption that Q is a definite germ cannot be omitted from the statement 



of Theorem 3.11 in view of Proposition y^^ v ) ■ Indeed, in the above proof, the fact 
that Q* is an germ is crucial. An example of non-uniqueness for an L^D germ 



which is not definite is given in Section 4.7 (Example 4.12) 



(ii) Let Q be a definite germ. In view of Proposition 2i.A(iv) and Definition 3.8 



it follows that Q- and Q* -entropy solutions coincide. Certainly, one can choose to 
work exclusively with maximal L^D germs (in which case Q* ^ G)- Yet there can 
be an advantage in using definite germs Q that are smaller than the corresponding 



maximal L^D extension Q* ; see in particular Sections^ pi |6.l[ and 6.2 



Remark 3.13. Let Q be a definite germ, and assume {u_,u^) Cz U x U is such 
that there exists a Q -entropy solution u of the corresponding Riemann problem 
(01,(0]). Then u is of the form (RPb-sol.) with {u\u'' ) £ Q* ■ In particular, 
u admits strong left and right traces on (0,cxd) x {0} which are equal to it' and 
, respectively. Indeed, the scaling invariance (Al) (see Remark 3.9(iii)) and the 



uniqueness statement in Theorem |3.11| imply that u is self- similar; by Definition 
\i.%\ i), it follows that u is monotone in the variable x/t on (— cx),0) and on (0,c»). 
Hence the traces u''^ := j^''^u exist, and by Remark 3.9(iv), the solution u is of the 
form (RPb-sol.) with [u^u") G g* . 

We want next to provide an equivalent definition of solution in which the trace 



condition of Definition 3.8 \i) is incorporated into the "global" entropy inequalities. 



As a preparational step, we state the following elementary lemma. 

Lemma 3.14. Let U be a closed interval in M and f E C(C/;M); for z,c G U, set 
q{z,c) — sign(z — c) (/(z) — /(c)) . Then for all a,b,c E U, 

\qia,c)-q{a,b)\<R{fi-y,c,b), (3.6) 

where i?(/(-); c, b) takes on one of the following forms: 

• -/?(/(•); c, 6) = 20sc(/(-); c, 5) , where Osc denotes the oscillation of f on the 
interval between c and b: 

Osc(/(-); c, b) :— max I min{c,6}<z<s<max{c,&}}; (3.7) 

• R{f{)\c,b) := 2u){\b — c\), where w is the modulus of continuity of f on the 
interval between c and b: 



to 



(h) := max | \f{z) — f{s)\ | min{c, b} < z < s < max{c, 6}, 1^ — s| < /i, } 



. provided fE BVi,,{U), R{f{-);c,b) 
f on the interval between c and b ); 



fc I /'(•*) I '^^ (this is the variation of 



Definition 3.15. Let Q be an L^D germ. A function u E L°°{R+ x M; U) is called 
a Q-entropy solution of (1.7), (1.2) if it is a weak solution of this problem and for 
all (c',c'') eU xU, 



[\u{t,x)~c{x)\^t + ^{x,u{t,x),c{x))^j,] - / \uo{x) - c{x)\^{Q,x) 



(3.8) 



+ / i?e((c',c'-))e(t,0) > 0, veel?([0,oo) xM), e>0. 
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where c{x) is the step function 
c{x) = c'n 



{x<Q} 



{x>0} 




the entropy flux q is defined by (13. 1| , and the remainder term Rg is given by 



Rg{{c\d-)):^2 inf (Osc(/'(.); c', 6') + Osc(r (•); c^ 6'^)), 
(f)',6'")eS 

with the oscillation function Osc defined in (13. 71) . 



Remark 3.16. The explicit requirement that u should be a weak solution is only 
needed to ensure that the Rankine-Hugoniot condition holds on S = (0, oo) x {0}. 
However, for many particular choices of Q , a bounded function u satisfying (3.8) is 
automatically a weak solution of (1.7), (1.2). 

Let us explain the choice of the penalization (remainder) term iJg ((c^, c'")) in 
Definition 3.15 which is inspired by an idea of Otto [S3]. The remainder Rg ((c', c'')) 
is chosen to satisfy the two foUowing properties: 

{c\c^)(.g^Rg{{c'.,c'-))=0, (3.9) 

and (because Q is an L^D germ) 

(c',c'^) e L/ X [/, (a', a'-) e ^ ^ q'{a\cT) - q\a\c^) < Rg{{c\c'-)). (3.10) 

Property ( 3.10[ ) follows from (3.6) in Lemma 3.14 via the bound 

for all {a\a''),{c\c'') eUxU, 

inf(,,,,.)ee|[g'(a',6')-'Z'^K,6")] - [<z'(a', c') - g'■(a^ c"^)] | < i?e ((c', c'')) , 

and the inequality q''{a\ b^)—q^{a^, ¥) > 0, which is valid for all (a', a^), (6', b^) e Q. 

According to Lemma 3.14 we can take the remainder term Rg in (3.8) under 
a few different forms, using, e.g., the moduli of continuity of the functions or 
their variation functions instead of the oscillation functions (3.7). If /'''' are globally 
Lipschitz continuous on U, the simplest choice is to take 

i?e((c',c'-)) 2 Wif'n'h^ inf - c'| + l^-- - | (6',6'') G G} 
= Cdist((c',c'^),g), 

where dist is the euclidean distance on and C is a sufhcienly large constant. With 
this choice, the properties (3.9) and (3.10) remain true. See [8^ for one application. 

The above Definition |3.15 is suitable for many generalizations, including the 
multi-D setting as in [21 [10] and the case of time-dependent families of germs as 
in [S]. However, for the precise model case (1.7), one can avoid the use of the 
remainder term Rg in (3.8), thanks to the following result. 



Proposition 3.17. The inequalities in (3.8) hold for all {c ,c^) lE U x U and for 

all nonnegative test functions ^, if and only if they hold (with zero remainder term) 
for the choices 



{c'' = c^ ^ c with c arbitrary and S,\x=o = 
U {(c', c*^) G G and ^ > is arbitrary} 



0} 



Proof. According to ( |3.9[ ), whenever {c\cT) £ G, the term i?g ((c', c''))^(i, 0) in 
(3.8) vanish es; cl early, it also vanishes if ^|a;=o = 0. Conversely, following the proof 
of Theorem 3.18 below, using the test functions in (3.11 ) and property (3.101 of 
the remainder TZg, it is easy to establish (3.8) for all choices of (c', c''), f > 0. □ 
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Using Proposition 3.17 with specific choices of G, from Definition |3.15| we recover 
the formulations of Baiti-Jenssen [M], Audusse-Perthame [12] (restricted to the 
model case (1.7)), and Biirger, Karlsen, and Towers [H]; see Section [1] for details. 
Namely, for test functions ^ which are zero on the interface {x = 0}, the Kruzhkov 
entropy inequalities with any entropy r/(z) — \z ~ c\, c £ R, are required; and 
for general ^, up-to-the-interface entropy inequalities are required only for a careful 
selection of "adapted" entropie^ The adapted entropies are x-dependent functions 
of the form 

T]{z,c) := \z - c{x)\, c(a::) = c' ll{a.<o} + c'' ll{^>o}, (c',c'')eC/. 

The notions of solution introduced in Definitions 13.81 and 13.151 bear the same 
name. Indeed, we have 



Theorem 3.18. For any L^D germ Q, Definitions 3.8 and 3.15 are equivalent. 



Proof. Consider a solution u of (1.7),(1.2| in the sense of Definition 3.8 and fix a 
pair (c', d") G U x U; consider now the function c{x) = ll{x<o} + ll{2;>o}. Take 
0<^g2?((0,oo)x]R) and consider the compactly supported in fi'^'' test functions 



{x^-h) + 
h 



(3.11) 



in the Kruzhkov entropy formulations of Definition 3.8 i). Since c{x) is constant in 
each of the domains fl''"^ , setting = + we have 

{\u - c{x)\ {£,h)t + qix,u,c{x)){^h)x} > 0. 



Sending /i > to zero, using Definition |3.8K i) for the existence of strong traces and 
Remarks |2.2l |2.3| for some of their properties, we deduce as in (3.4) that 



(3.12) 



{\u{t,x) - c{x)\^t + q{x,u{t,x),c{x))(x} 

{Q'{i-f'v'u){t),V'c')^Q^{{YV^u){t),V^-c^))am>0. 



Set w'''^(t) := (7'''V''''u)(t) and z'''^ c''''. By Definition [3!8[u), 

{v^(t),v''{t)) eVg*, fora.e. i>0. 



Hence, by Remark 3.7 V(?«',w'') £ VQ we obtain Q\v^,w^) > Q'"(ti'', w*"), and so 
Q\v',z') - Q'■(^;^ z'') = Q\v',w') - Q^(^;^ w'') + Q\v\z') - Q'\v\ z'') 

-Q\v\w') + Q'^{v\w^) 
> [Q\v\z')-Q^{v\z^-)\ - [Q\v'M)-Q'-{v\wn]- 

coming back to the definitions of V'^'^ and Q''^ From ( |3.10[ ) we then easily deduce 
that 

Q\v\z')-Q^{v\z^)>-Rg{{c\d))- 

therefore (3.12) yields (3.8), at least for test functions ^ that vanish on {t — 0}. 

We deduce ( 3.8 1 for all < ^ G I?( [0, oo) , M) , by truncating ^ in a neighbourhood 
of {t = 0} and using Definition |3.8[iii) . Thanks to Remark [3. 9[ii), u is also a weak 



solution solution of (1.71, (1.2); thus it is a solution in the sense of Definition 3.15 



These formulations are very similar to the one introduced by Carrillo in 1221 for homogeneous 
Dirichlet boundary value problems for degenerate parabolic problems including conservation laws. 
In our context, we can use the standard Kruzhkov entropies \z—c\ in the place of the semi-Kruzhkov 
entropies [z — c)^ used by Carrillo. 
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Reciprocally, take a solution u of ( |1.7[ ), (1.2) in the sense of Definition 3.15 
First, it is clear that Definition [S^Sf^i) holds. Therefore it follows from (2.2 1 that 
also Definition [SjSjiii) is fulfilled. Similarly, it follows from Theorem 2.1 that strong 
left and right traces 7''''^''''^ on E = (0, 00) x {0} exist. It remains to show that 
{{j^V^u){t), {YVu){t)) e Vg* pointwise on S. To this end, with ^ e X'([0,oo) x 



Hi), C > 0, and ft. > 0, take i}}^ = |^^^^ ^ as test function in (3.8 1 and pass to the 
limit as /i — > 0. Using Remarks |2.2| and |2.3[ we deduce 

{Q\{^'v'u){t)yc')-Q^{{YV^u){t),V^c^))am>Q. 



Now we take (c',c'") G Q and utilize that i?g ((c', c'')) = 0, cf. property (3.9). Since 
^(•,0) e 2?([0,oo)), C > 0, is arbitrary, we get 

Q\{-i^V^u){t)yd) >Q''{{-i'^Vu){t),Vc''), fora.e. t>0. 



for any (c', c'') € Q. In addition, recall that, u is a weak solution of ( 1.7); by Re mark 
the Rankine-Hugoniot relation (7'(7'F'm) = g"^ {'-^'^V'^u) holds. By Remark 



2.3 



3.7 



we infer that 



This justifies Definition |3.8[ii) and concludes the proof. 



{{-i'v'u){t),{^^v'u){t)) evg* 



□ 



3.3. Comparison and continuous dependence results and L°° estimates. 

In this section we provide a comparison result and some estimates for tj-entropy 
solutions. Moreover, we introduce a "distance" between two germs, which is used 
to state a continuous dependence result with respect to the choice of the germ Q. 
This result prepares the ground for a study of parametrized families of germs in [3] . 

Theorem 3.19. Assume that Q is a definite germ. Let u and u be two Q-entropy 



solutions of problem (1.7), (1.2) corresponding to initial datauo anduo, respectively, 
such that {uq — uq)^ G L"'"(M). Then 



{u ~ u) (t) < / (uq — iio) , for a.e. t > 0. 



In particular, uq < uq ^ u < u. 



Proof. The proof is the same as the one used to conclude Theorem |3.11| The 
doubling-of- variables arguments is used to derive the Kato inequality with (u — w)"*" 
instead of |w — u|, and we then derive the an alog ue of inequality ( |1.11[ ) for the 

these fluxes are expressed by 



2.2 



associated entropy fluxes q'^ . As in Remark 
means of the functions 

Q'i' : (2, k) ^ sign±(z - k)ig'-'^{z) - g'^^k)). 

Take {v',v''), {v\v^) e VG*; it suffices to show that 

By Remark 3.7 we have (3.3) and, moreover, g'(w') 



(3.13) 

g^ivn, gKv') = g^{vn- 



Consider for example the case g^{v^) — > 0, in which case (3.3) is equivalent 

to the inequality sign(w' — w') > sign(i;'" ~ v^). This implies the following two 
inequalities: sign^(t;' —■()')> sign*(u'' — ■0'"), which in turn imply (3.13). The other 
cases are similar. □ 

As a corollary, we have the following maximum principle: Let u be a C/-entropy 
solution of ( 1.7) with an initial function uq satisfying for some (c', c""), (C', C) S Q* 



the following lower and upper bounds: 
< uq{x) < & for a.e. a; G IS- 



C'' < un{x) < C for a.e. x S K+. 
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Then for a.e. t > the same inequahties are satisfied by u{t,-). The proof is 
immediate from Proposition 3.10[ ii) and Theorem 3.19 

In general, \\u{t, can be greater than |1uo|Il°°J ^-S- the examples in pU] . 

But we have 

Proposition 3.20. Let Q be a definite germ. Let ni^ < M_, m+ < M+ be 

real numbers such that there exist Q -entropy solutions of the Riemann problems 
(1.7), (1.8) with data = (m_,m_|_) and with data = (Af_,M+). 

By Remark 3.13 strong left and right traces of these solutions on {x — 0} exist; 
denote these traces by c'''',C'''", respectively. If m± < uo < M± for a.e. in and 



u is the Q-entropy solution of (1.7) with the initial function uq, then 
min{m_, c'} ll{:r<o} + min{TO+, c''} ll{j;>o} 



< u{t,x) < 
max{M_, C'} ll{,<o} + max{A/+, C"} 11{, 



{x>0}i 



Proof. The claim follows from the comparison principle in Theorem 3.19 and from 
the monotonicity property of the solutions of the Riemann problems in the domains 
{a; < 0}, {a; > 0}, cf. Remark [313} □ 



Finally, we state a simple result regarding the continuous dependence of G- 
entropy solutions on the choice of the L^D germ Q. For a pair of fixed functions 
/'''" G C(J7;M) and the associated Kruzhkov entropy fluxes consider two L^D 
germs Q, Q. Define the "distance" function 



:-max<^0, sup {q^- {b\b'') ~ q\b\b')) 



(3.14) 



Observe that p{G', G) = implies G C G* and G C G* ■ So if the distance p between 
two maximal L^D germs is zero, then the two germs coincide. 

Proposition 3.21. With the notation above, let u, ii be G- and G-entropy solutions 



of (1.71, respectively, with initial data uq,uq, respectively. Then 



\u — u\{t)< / \uq ~ uq\+ t p{G,G)-, fora.e.t>0. 



(3.15) 



The proof of Proposition |3.21| represents an obvious modification of the last 
argument of the proof of Theorem 3.11 (namely, the term (3.5 1 is controlled by the 



last term on right-hand side of (3.151 when < ^ < 1 



3.4. Complete germs and tj-entropy process solutions. In this section we 
will discuss briefly the question of existence of t^-entropy solutions for a general 
definite germ G, and introduce the weaker notion of tj-entropy process solutions as 
a tool to establish convergence of approximate solutions equipped with mere L°° 
bounds. For some particular cases, existence will be shown in detail in Sections |4] 
and 6.1 while more general results will be given in Section [6] 



Definition 3.22. A germ G is said to be complete if for all (w_,u+) dU xU, the 



Riemann problem (1.7), (1.8) admits a weak solution of the form (RPb -sol .) with 



G G (here and in the sequel, we adopt the convention of Remark 1.2) 
Observe that by Proposition 



3.4 



i) and Remark 3.9 ^iv), if G is an L^D germ, 
then the aforementioned weak solution of problem (1.7), (1.8) is also the unique 
C/-entropy solution of the problem. 

Proposition 3.23. Let G be a definite germ. If there exists a G-entropy solution 
of (1.71,(1.2) for any Uq G i°°(M; U), then G* must be complete. 
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Proof. Consider general Riemann initial data of the form (1.8). By assumption 



and by Remark 3.13 there exists exactly one C/-entropy solution u to (1.7), (1.8); 



moreover, u is of the form (RPb-soL), and the corresponding pair of states (m', u^) 
belongs to G* . Since u is necessarily a weak solution, Q* is complete. □ 

The assumption that G* is complete is expected to yield existence of tj-entropy 
solutions. Indeed, starting from an appropriate Riemann solver (a classical solver 
for X ^ 0, and the solver TZS^ at a; = 0), it is possible to construct numerical 
approximations of C/-entropy solutions, e.g., by the Godunov finite volume scheme 
or by the front tracking scheme; see Section |6.3[ Unfortunately, currently we can 
deduce existence only in the presence of good pre-compactness properties (see, e.g., 
[35] and [21] for uniform BV estimates for front tracking and adapted Engquist- 
Osher schemes, respectively). 

In passing, we mention that "completeness" makes available uniform L°° bounds. 
If Q is definite and G* is complete. Proposition 3.20| yields an L°° bound on any 
^-entropy solution SfuQ of the form 

||(5f uo)(-)||i=e(i{) ^ Const (||uo||l~,^^) , uniformly in 

If similar L°° estimates are available for sequences produced by an approximation 
procedure, and if the existence of the tj-entropy solution is already known by other 



means, then Definition 3.24 and Theorem 3.28 below provide a "propagation of 
compactness" approach to proving convergence of such sequences. 

To this end, following [TH IHS] [331 IM] (see also references cited therein and the 
previous works of Tartar and DiPerna) , we first extend the notion of solution given 



in Definition 3.15 to account for so-called process solutions (higher dimensional L° 



objects related to the distribution function of the Young measure). 

Definition 3.24. Let G he an L^D germ. A function fie L°°(R+xRx (0, 1); U) 



is called a G -entropy process solution of (1.7)Jl.2| if the following conditions hold: 



1 (weak process formulation) . For all ^ G 2?([0, oo) x M), 



{pi(t, a;,Q;)^t + f(x,^(t, a;,Q;))Ca;} - / uo^(0,a;)=0 



(3.16) 



2 (penalized entropy process inequalities). For all pairs {c}^c^) £ U x U and for 
all ^ eV{[0,oo) xE), ^>0, 



{\li{t,x,a) - c{x)\^t + ci{x;tJ.{t,x,a),c{x))^^} 

\uo ~ c{x)\aO, x)+ f Rg ((c', c'-))e(t, 0) > 0, 



(3.17) 



where c(x) = c^ ^{x<o} -^c^^ixyo}, o,nd the remainder term _Rg ((c', c*")) has the 
same meaning as in Definition \'i.lb\ 



As for Definition 3.15 alternative (equivalent) forms of the remainder term can 



be chosen, provided the two properties (3.9) and (3.10) are fulfilled 



Remark 3.25 (cf. Proposition 3.17). A G-entropy process solution is equivalently 
characterized by the following three requirements: 



(a) The weak process formulation (3.16) holds; 

(b) For all ^ G V{[0, oo) x M), ^ > 0, such that ^ = on the in terface {x = 0}, 
for all pairs (c, c), c £ U, the entropy inequalities ( 3.17[ ) hold (with zero 
remainder term); 

(c) For all ^ G X'([0 oo) x M), ^ > 0, for all pairs (c',c'') G G, the entropy 
inequalities (3.17) hold (with zero remainder term). 



ON SCALAR CONSERVATION LAWS WITH DISCONTINUOUS FLUX 



19 



Indeed, it is evident that Definition 3.24 implies the three properties above. 
Clearly, in the cases (6) and (c) above, the remainder in (3.17) vanishes, cf. (3.9). 
For the proof of the converse implication, as in the proof of Theorem |3. 18 we utilize 
the test functions given in (3.11) and the remainder term property (3.10). 

Remark 3.26. Fix a pair (c', c'') € UxU, and let fi be a Q- entropy process solution. 
We want to consider the functions 



I 



Lr 



{t,x)i-^ I \fi{t,x,a) - c'''''\ da, j''*" : (i, x) / x, a), c'''') da. 



J,r\ 



The inequalities in (3.17) imply that the vector field 

[t,x)^ {I\t,x),j\t,x)) 

is an L°° divergence-measure field on M-|_ x (— oo, 0) (locally in t). According to the 
results in Chen and Frid |24j , such a vector field admits a weak normal trace on the 
boundary S = K_|_ x {0}. This means that j' admits a weak trace from the left on 
E. Similarly, we conclude that admits a weak trace from the right on S. 



3.26 



In what follows, the weak trace operators t i— > 7^'' J •) from Remark 
will be denoted by jlf . 

Lemma 3.27. Let ji be a Q -entropy process solution. Then for all {v\ v^) G VQ , 



iL Q\v'pL(t,;a),v')da>^l, Q' (y>(t, .,«),«'■) da 



(3.18) 



for a.e. t > 0. 

Proof. The proof is very similar to the proof of the second part of Theorem |3.18[ 
We take c'''" such that 

we take M)_^ as test function in 



= V^-''d'\ For i e ^([O, 00) X M), C(-, 0) > 0, and ft > 0, 



the limit as ft to finally arrive at (3.18 



3.17[ ), use the Fubini theorem, and pass to 

□ 



Notice that, unlike Remark 3.12| [ii), it is not clear whether Q- and Q*- entropy 
process solutions coincide for a definite germ Q . Therefore we require that Q = Q* 
(i.e., we require that Q is a, maximal L^D germ) while working with tj-entropy 
process solutions. 

Theorem 3.28. Let Q be a maximal L^D germ. Suppose Uq is such that there exists 
an Q-entropy solution u of (1.7| with initial data (1.2). Then there exists a unique 
Q-entropy process solution /i with initial data uq, and ^(a) — u for a.e. a G (0, 1). 



Proof. Bein ^ a ^-entropy solution, u fulfills Definition |3.8[ ii). Since G = G* 
Proposition 



3.4 



;iii), the pair {v\t),v'- {t)) := {{-f^V^u)it), {YVu){t)) belongs to 
VG* for a.e. t > 0. Hence, inequality (3.18) at a point t > holds in particular 
with the choice w'"'" — v'''^(t). Further, notice that since the traces of V'^'^u are 
strong, the weak traces jlf Q'"^ {V'^'^ fi{a), V'"^u) da on S exist and are equal to 



Jo i9''''(V'''Xa),7'''V''''u) da. Therefore 



{V^ti{t,-,a),V^u{t,-)) da 



7; / Q''{V-tiit,;a),V''u{t,-))da]^{t,0)>0 



(3.19) 



for all ^(-,0) e T>{[0, 00)), £,{■, 0) > 0. Now we repeat the arguments in the proof of 
Theorem 3.11 which leads to the analogue of inequality (3.4); after that, in view 
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of (3.19 1, we deduce 



{\u - fi{a)\£,t + q{x,u,^{a))£,^} > 0, 



because u and fj, are solutions with the same data. As in [50| (for Lipschitz fluxes 
/'''') or in [iniET] (for merely continuous /'■''), we work out finally that 



^(•, •, a) = u{-, •) a.e. on x M, for any a G (0, 1). 



□ 



We have the following open problem: Is it true that for a definite germ Q there is 
at most one Q -entropy process solution fi with given initial data Uq (in which case 
fi is independent of a) ? The answer is known to be positive in some situations. 
Firstly, according to [331 131], this is the case for the Kruzhkov entropy process 
solutions of ut + f{u)x = 0, which can be eq uivalently defined by Definition 3.24 
starting from the germ QvKr (see Section 4.2). Secondly, if /' is non-decreasing and 



is non-increasing, then every L^D germ ^ is a singleton {(«',«'')} (see Section 
4.5); in this case, from (3.181 we readily derive that V^'^ ^{a) = u'''' for a.e. a, and 



thus uniqueness of a tj-entropy solution follows. 



3.5. Conclusions. Let us review the results of Section|3] Admissibility of entropy 
solutions to the model equation (1.7) is most naturally defined starting from a germ 
Q possessing a unique maximal LrD extension Q* that is complete. In this case, 
uniqueness, L^-contractivity, and comparison properties hold for solutions in the 
sense of (the equivalent) Definitions 3.8 and 3.15 and we can hope for the existence 
of solutions for general data. Furthermore, uniform L°° estimates of solutions can be 
obtained (see, e.g.. Proposition 3.20) and measure-valued solutions (see Definition 



3.24) can be considered, but for the latter we need to know already the existence 



of solutions in the sense of Definitions |3.8[ |3.15| and also the additional assumption 
Q — Q* must hold; under these conditions we are able to conclude that Definition 
|3.24| is in fact equivalent to Definitions 3.8 3.15 In Sections [4j [5] we illustrate the 
theory of germs by a number of examples. 

Remark 3.29. Each maximal L^D germ Q gives rise to an -contractive semi- 
group for (1.7). Conversely, by the analysis in Section \l.2\ any semigroup defined 
on the whole space L°°{R; U) that satisfy (Al) and (A2) corresponds to a complete 
maximal L^D germ. If we drop the semigroup property and the scaling invariance 
(Al), then different L^-dissipative solvers can be constructed, e.g., starting from 
a family {G{t)}^^^ of maximal L^D germs; we refer to Colombo and Goatin [26] 
for an example of such a solver, and to |5] for an analysis of the Colombo-Goatin 
solver in terms of admissibility germs. Notice that in order to have existence we 
must assume that Q{-) is measurable in an appropriate sense. We will develop an 
approach to measurability for time- dependent families of germs in [9]. 



4. Examples and analysis of known admissibility criteria 

In this section we first discuss the completeness of germs, a "closure" operation 
on germs, and relations between completeness and maximality of L^D germs. We 



continue by reviewing a number of known admissibility criteria for problem (1.7) 



and make explicit the underlying germs. In some cases, using the results of Section 
|3] we refine the known uniqueness theory. Some known and some new existence 
results for these criteria are contained in Sections |6.2| |6.3l A similar study of the 



important "vanishing viscosity" germ is postponed to Sections [5] and |6.1[ 
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4.1. More about complete germs and closed germs. For a given G U, we 
denote by 6''(-, u') the graph that contains all points /'(m_)) G [/ x M such that 
there exists a Kruzhkov solution to the Riemann problem joining the state u_ at 
X < 0, t = to the state u' at x — 0~ , t > (i.e., this solution contains only waves 
of nonpositive speed). Similarly, for u*" G U, 9'^{-,u^) is the graph of all points 
/''(«+)) such that there exists a Kruzhkov solution to the Riemann problem 
joining u'' at a; = 0+, i > 0, to at a; > 0, t = (i.e., this solution contains only 
waves of nonnegative speed) 

Notice that in concrete situations, the completeness of a given germ can be 
checked as follows. 

Remark 4.1. A germ Q is complete if and only if 

IJ Dom6i'(-,?i') X Domr(-,it'') = C/ X [/. 



According to the convention of Remark 1.2 the pair of traces (7'u,7''u) of a 
solution of the form (RPb.-sol) does not necessarily coincide with [u^^u^). The 
reason is that zero-speed shocks may be, at least in principle, a part of the wave 
fans joining m_ to u' and to . In order to cope with such situations, we 
introduce the following definition. 

Definition 4.2. A left (resp., right) contact shock is a pair (u_, u') E U xU (resp., 

a; < 

,' ' is a Kruzhkov entropy solution 

u , X > 

{u^ X < 
' ^ is a Kruzhkov 

u+, X > 0, 

entropy solution of Ut + f^iu)^ — in x Mj. 
In other words, 

{u^,u') is a left contact shock ■^=> w_ G Dom 0\-,u^) and f\u-) — f{u^); 

an analogous statement is true for the right contact shocks. 

Clearly, any pair (c, c) G [/ x U is both a left and a right contact shock. 

Definition 4.3. A germ Q is said to be closed, if Q is a closed subset of U x U 
and, moreover, for all pairs (u^u^) E Q, G also contains all pairs (m_,w+) such 
that {u^,u^) is a left contact shock, {u^,u-^) is a right contact shock. The smallest 
closed extension of G is called the closure of G and is denoted by G- 

It is clear that for any solution u of the form (RP b.-so l) such that (u^u^) G G, 
we have {'~f''u,j^u) G G (cf. the convention of Remark 1.2). 



The above definition is consistent in the sense that the closure C/ of a germ G 
is indeed a germ. Indeed, the Rankinc-Hugoniot condition for {u-,u^) G G holds 
because it holds for the corresponding pair {u^u"^) G G- 

Let us list some properties of the closure operation and of closed germs (to be 
utilized in Section [s]). 

Proposition 4.4. 

(i) G is an germ ^ G is an L^D germ. Furthermore, G G* , {G) = G* ■ 

(ii) A maximal L^D germ G is closed. 

(Hi) Any maximal L^D extension of G contains G- In particular, G is a definite 
germ if and only if G is a definite germ. 

(iv) If G is a definite germ, then G and G entropy solutions coincide. 
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Proof, (i) Fix {u\u^), {u^vT) S Q, and let {ii-.v}) be left contact shocks 

and (u'', u+), {u^ , u+) be right contact shocks. The Kruzhkov admissibility of these 
shocks implies that 

Moreover, we have q^{u\v!-) > q^{u^,u^) because G is an L^D germ. Therefore we 
infer q^{u^,u^) > q^ {u^,u^). From the definition, it is not difficult to see that Q 
is obtained by the following two operations: 

- firstly, taking the topological closure in M^, and 

- secondly, adjoining (possibly trivial) left and right contact shocks to the 
elements of G- 

The result of the second operation is topologically closed. Because each of the two 
operations preserves inequality , by the argument above and by the continuity 



of (?'''', G is an L^D germ. By Proposition 
G is contained within G* 
Finally, by Definition 



3.4 



ii), G being an L^D extension of G, 



3.2 



G C G implies (G)* ^ {G) ■ The reciprocal inclusion 
ii) , each pair {u\u^') G G* belongs to an 



3.4 



is also true. Indeed, by Proposition 
L^D extension G' of G- By Definition |4.3[ G C G' implies G C G'- By our first 

ii), we infer 



3.4 



claim, G' is an L^D extension of G, using once more Proposition 
that {u^u"-) € Wc (G)*. Thus G* = (G)*. 

(ii) By Definition |4.3| and by (i), one has G Gl G Gl G* ■ Since t/ = t/* by Proposition 
|3.4[ iii), it follows that G coincides with its closure. 

(iii) li G '^G' and G' is a maximal L^D germ, then G C G'] besides, G' — G' by (ii). 

(iv) This follows from (i),(iii) and Remark 3.12[ ii). □ 

Let us list separately the properties related to complete germs. 
Proposition 4.5. 

(i) If G is a complete L^D germ, then G is a maximal L^D germ. 

(ii) If G is an germ such that G is complete, then G is definite and G* = G . 

Notice that in case (ii), ^ is a definite L^D germ and G* is complete; as it is 
pointed out in Section |3.5[ such germs are expected to lead to a well-posedness 
theory for -entropy solutions. 

Proof, (i) Let G' be a nontrivial maximal L^D extension of G, and pick a pair 
(u_,u+) e G'\G- Now the Riemann problem (1.7), (1.8) possesses two solutions of 
the form (RPb.-sol), namely 

- a solution u — TZS^ {u-, u+), which exists because G is complete; 

- the elementary solution u = u- ll{2;<o} + u+ ll{a;>o}- 

iv) and because {"f^u,^^) G G C 



3.9 



Both are C/'-entropy solutions, by Remark 

G' = {G')* and (7'u,7''-u) — (m_,u+) e G' = {G')* ■ In addition, u and u do 
not coincide a.e. because their traces are different: (7'M, 7''u) = (u_,u+) ^ G ^ 

(7*^,7'' u)^ Since G' is maximal L^D, this contradicts the uniqueness result of 

Thi s con tradiction proves that G is itself maximal L^D. 

i), is a complete L^D germ. By (i), it is a maximal L^D 



Theorem 



3.11 



4.4 



(ii) By Proposition 

germ and th us it is def inite . By Proposition 4.4 |iii), also G is definite. Finally, by 
Propositions 4.41) and 3.4111), it follows that G* ^ (G)* ^G. □ 



Finally, let us point out a situation where the existence of a specific pair in an 
L^D germ G immediately excludes some other pairs from G- 

Proposition 4.6. Let G be an L^D germ, and fix a pair {u'',u^) e G- Consider 
e Dom6l'(-,u') and uj\- G Dom 9^{-,u^). Then (u_,u+) ^ G, except in the case 
with f\u^) = f-{u+) = f\u') = f^{u-). 
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Proof. By Proposition 4.4 'i), Q is an L^D germ. Since G G, in case 

E Q we must have f{u^) — and q^{u^,v!') > {u^,u^). But 

the monotonicity of 6''^''(-, u'^'') (which is easy to get from the description of the 
Kruzhkov solutions of the Riemann problem, see, e.g., [40l|39]) yields 

q\u^.u') = sign(u_ - u'){f{u^) - f{u')) = - \f{u^) - f{u')\ , 

q'-iu+y) = Sign(u+ - u^Wiu^) - fiu^) - |.r(u+) - rK)| . 

These conditions are not compatible, unless f{u^) ~ f^{u^) = f^{u^) = /' (u''). In 
addition, if the latter condition holds, then (resp., is a left (resp., 

right) contact shock; in which case E Q. □ 

4.2. The case /' = /''; the Volpert-Kruzhkov germ. We will now illustrate 
the abstract "discontinuous flux" theory on a well-known example, namely the case 



without a flux discontinuity at {x = 0}, so /' = in (1.7); let us write / instead 
of /'-'■ and take U = R. 



In this case, the set of weak solutions of (1.7) includes all constant in M+ x M 
solutions. Let us postulate that all constants are admissible elementary solutions 
of ut + f{u)x = 0; in other words, consider the "Volpert-Kruzhkov" germ 

GvKr :-{(c,c)|cGM}. 

One checks immediately that G is an germ. According to the definition, the 
dual germ Gyxr determined by the family of inequalities: 

{u',u-)eG*vKr^q{^',c)>q{u-,c), VceM (4.1) 

(here q{u,c) — sign(u — c){f{u) — f(c)) ). In deed, the Rankine-Hugoniot condition 



/(m') = f{u^), which enters Definition 3.2 follows by taking ±c large enough in 
(4.1). Following Volpert [HHl: we check that Gyxr an L^D germ: 

q{u'- ,v}) > q{vJ' .v!-) = q{u'-,u'') > q{vJ' .vT) = q{vJ',vJ'), 



for all {u ,u^), {u ,vT) € Gyxr- Proposition 3.4 V), GvKr is definite. Therefore 
uniqueness of C/ i/ifr-entropy solutions holds. 

Moreover, it is well known (see, e.g., and [JOl [SH]) that there exists a 



Kruzhkov entropy solution of ut -I- f{u)x = with the Riemann data (1.8 1, and 



that this solution is of the form (RPb.-sol) with {v} ,u^) satisfying (4.1). Therefore 
the germ G*jij. is also complete, which raises the hope for existence of yifr-entropy 
solutions for all bounded data. Not surprisingly, of course, it easily follows that 



Definition 3.15 with G = GvKr is equivalent to the definition of Kruzhkov entropy 
solutions. Hence, the existence of i/ifr-entropy solutions is a consequence of [50) . 

In conclusion, the germ GvKr leads to a complete well-posedness theory for 
solutions in the sense of Definitions |3.8| and |3.15[ which turns out to be precisely 
the theory of Kruzhkov entropy solutions. Another description of the associated 



maximal L^D germ GvKr is given in Section 4.3 below 



Remark 4.7. As a final remark, this example is somewhat misleading. Indeed, we 
have just investigated the question of admissibility of discontinuities of solutions of 
Ut + f{u)x = Q on the line {x = 0}; and we took for granted that the admissibility 
of u on {x < 0} and {x > 0} is understood in the sense of Kruzhkov. The focus on 
the line {x = 0} is "artificial" in this example. The fact that we finally arrive at 
notion of Kruzhkov entropy solutions on the whole o/M_|- x M is o/ course expected; 
moreover, our existence analysis relies on the knowledge of the Riemann solver 
associated to Kruzhkov entropy solutions! But we believe that this example and its 
sequel in Section \A.2>\ illustrate well the general approach and some techniques related 
to the use of admissibility germs. Also see Section \A.9\ and iZemarfc |4.8| /or examples 
of non-Kruzhkov L^-dissipative solvers for Ut + f{u)x — 0. 
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4.3. The case /' = /''; the Gelfand germ. We will continue with the above 
example but apply a different approach, which will be extensively investigated in 
Section [5] in the case of a discontinuous flux. 

Recall that existence of Kruzhkov entropy solutions for 

ut + f{u)x = 0, u|f=o = ua (4.2) 

was first obtained in [50^ through the limit, as £ ^ 0, of the vanishing viscosity 
approximations 

Ut + f{u)x = euxx- (4.3) 



Indeed, it is well known that there exists a unique solution u'^ to (4.3 ) for all bounded 
data, that the solutions satisfy an analogue of the L^-dissipativity property (1.5) 
for each e > 0, and that for any uq € L°°, (u'^)e>o is relatively compact in Ll^^. 



Using diagonal, density, and comparison arguments (see Section 6.1 below for 



details), at the limit e — > wc arrive at an L-'^-dissipative solver for (4.2 ), defined for 
all Uq € L°°. By Remark |3.29[ this implicitly defined "vanishing viscosity solver" is 
characterized by some complete maximal L^D germ that we now want to determine. 



Its dual germ Q* consists of all elementary solutions (1.10 1 that are trajectories of 
this solver. According to Gelfand j36] (cf. [40l[39]), for all (u'^u^) satisfying 

f(it') — f(u^) and either u' = , or 

(4.4) 

sign(M'" — v}){f{c) — f{u')) > for all c between u' and , 



there exist standing-wave profiles : M — )■ M such that lim^^_oo = W 



lim^^^oo W{i) = vT , and W{^) solves (4.3) for all £ > 0. Thus solutions of the form 
( |1.10[ ) with {u'-.u') as in (|4.4|) are explicit limits as £ -> 0, of solutions of (|4.3|. 



Therefore these solutions are trajectories of the vanishing viscosity solver, and we 



can define the "Gelfand" germ Qg as the set of all pairs (u^u^) satisfying (4.4) 
Gg is L^D because the vanishing viscosity solver is dissipative by construction. 

Let us look at the closure Qg of Gg- One checks that Qg D G' , where G' is the set 
of all pairs {u\u''') such that ( |4.4[ ) is satisfied with the inequality "> 0" instead of 
"> 0". But G' is known to be complete (see [331101131]), therefore Gg is a complete 
germ. By Proposition |4.5[ ii), Gg is a definite germ with complete closure; more 
exactly, Gg — Gg — G' , because G' turns out to be an L^D germ. Therefore we can 
identify the "vanishing viscosity" germ with the unique maximal L^D extension of 
Gg-i i-6- with its closure Gg- 

Now Gg- and t/ y/fr-entropy solutions must coincide, because GvKr C Gg (see 
their definitions) and GvKr was shown to be definite; in addition, we infer that 
G*vxr coincides with Gg- In conclusion, the germs GvKr a-nd Gg have different 
motivation, but they both lead to the classical Kruzhkov entropy solutions. 

Remark 4.8. In closing, we mention that non-classical -dissipative solvers exist 
for (4.2). For example, for the Hopf-Burgers flux f{u) = , Gg 7^ K x M. Take 



any pair (u',u'') ^ Gg, which satisfies the Rankine-Hugoniot condition, and then 
set G" = {{u^u'' )}. Take an arbitrary maximal D -extension G'" of G" ; the 
associated G'" -entropy solutions are, in general, different from the Kruzhkov entropy 
solutions! Notice that since (u^u^) ^ Gg = OvKn some of the constant elementary 
solutions in GvKr will not be trajectories for a non-classical solver. Finally, we 
recall that physically motivated — but not -dissipative — non-classical solvers for, 
say, the nonconvex conservation law Ut + (u^)x = have been extensively studied 
in the recent years, cf. the book of LeFloch I54j . 

4.4. Increasing surjective fluxes Baiti and Jenssen [H] have proposed an 



entropy formulation for (1.7) with strictly increasing and uniformly Lipschitz con- 
tinuous fluxes /' , f^' (strictly decreasing functions can be treated similarly) . 
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In what follows, we assume that there is a set C M such that /'''' : U W 
are non-decreasing and onto. The corresponding germ Grh consists of all pairs 



(u ,u^) €U X U that satisfy the Rankine-Hugoniot condition (1.9). 

Due to the monotonicity of /'^'', ( 1.11 ) holds for all (u^u^), {u' ,u^) e Qrh- Thus 



Grh is an L^D germ; it is a maximal L^D germ because Qrh admits no nontrivial 
extension. Finally, Grh is complete. Indeed, taking d U xU^ there exists 

& U such that [u-.u^) G Qrh\ take any e Furthermore, the 

graph 6^{-,u^) coincides with the graph of /''. Therefore G Dom 6'''(-, w''), and 
the Riemann problem ( 1.7 1, ( 1.8 1 admits a weak solution of the form (RPb-sol.) with 
the intermediate states G Grh- Definition 3.15 of C/^j/f-entropy solutions 

then turns into the definition of an admissible solution in the sense of Baiti and 



Jenssen JJj, adapted to the model equation (1.7). 

In passing, notice that a simpler equivalent choice of the admissibility germ is 
possible. Indeed, we can take Gii ~ 0- It possesses a unique max imal L ^D extension, 
namely (G^)* = Grh- Thus Gtii is a definite germ; by Remark 3.12 the Grh- and 
C/0-entropy solutions coincide. 



The entropy dissipation property (3.8) of an fjg-entropy solution u imposes no 



restriction on the jump of u across the flux discontinuity, so the only restriction 
is the Rankine-Hugoniot condition: 7'/'(u) =l^f^{u) (which is equivalent to 
jWu) e VGrh = V{Gii)* of Definition 

4.5. Other cases with monotone fluxes 



3.8 



ii), since /'''' are monotone). 

We have two different situations. 
If we assume that on [/ C M, /' is non-decreasing and is non-increasing, 
then the germ Grh of Section 4.4 is still a maximal L^D germ, but it is not 
complete. Indeed, for all m'"'" £ U the graphs 9^'^{-,u^-'^) are constant; they 
are the sections of the graphs of f'^ at the levels /'■''(ii''''), respectively. 
Therefore, unless the Riemann data (u-, m+) belong to Grh, there exist no 
solution to ([r7|,([r8| under the form (RPb-sol.) with {vJ-.u') € Grh- 
Consider the case where /' is non-decreasing and is non-increasing, on 
U CR. If fiU) n f{U) = 0, then the Rankine-Hugoniot condition 



never holds, i.e., there exist no weak solution of (1.7 1. 
Let us assume in addition that the ranges of /' and are not disjoint. We claim 
that in this case, for any pair (A, B) such that the Rankine-Hugoniot condition 
f''{A) = f^{B) holds, the germ G{a.b) — {(^7 B)} is a definite complete germ. The 
corresponding maximal L^D extension of G{a.b) is its closure G{a,b)- Indeed, it is 
evident that G(a,b) is an L^D germ. If t/ is a maximal L^D extension of G{a,b): 
then {A,B) S G- Notice that the graph 6'{-,A) coincides with the graph of 
and the graph 9^{-,B) coincides with the graph of /''. In particular, both graphs 
have U for their domain. It follows by Propositions 4.6 and 4.4 iu) that for all 



(w_,u+) gU xU, either (u_,u+) ^ G = G 01 (m_,u+) G G(a. b )- T hus G = G(a,b)- 
Finally, for all (u_,u+) E U x U, the Riemann problem (1.71,(1.8) admits a weak 
solution of the form (RPb-sol.) with (u^u^) 
Notice that the germs G{a.b), G(A'.B') 



{A,B). 

described above have different maximal 
L^D extensions if and only if f{A) = f{B) ^ f{A') = fiB'). Hence, whenever 
f^{U) n f^{U) contains more then one point, there exist infinitely many two by two 
non- equivalent definitions of ^-entropy solutions, each definition being based on a 
complete definite admissibility germ. 

4.6. The Audusse-Perthame adapted entropies. Now we consider the case of 
two fluxes /'■'' such that for some G f'''^ ■ {—oo,ul^^] !-> are bijective 
decreasing functions and : [Uq''',oo) i-^ M+ arc bijective increasing functions, 
cf. Figure [I] for an example. This is the context of Audusse and Perthame in [T^], 
reduced to our model case (1.7 1. The work [T^ applies in more generality; it was 
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recently revisited by Panov [7D], who showed that a suitable change of variables 
reduces the Audusse and Perthame definition to the definition of Kruzhkov [50] . 
Here we do not take advantage of Panov's observation. 



/ 














u 



Figure 1 . Typical left and right fluxes covered by the framework 
of Audusse and Perthame; Bold dots represent Wq'''. 



The notion of entropy solutions in [T^] is based on "partially adapted entropies" . 
These are, for our model case, entropies built from a step function c{x) associated 
with a pair (c',c'') G Qap'- z i->- rjc(x){z) ^ \z — c{x)\. Here Qap is the "Audusse- 
Perthame" germ given by 

{u\u^) e Gap ^ f in') - TK), sign(ii' - ui) = signlu"- - <). (4.5) 

Definition 3.15| (under the form of Proposition 3.17) of CJ^p-cntropy solutions 
reduces to the definition of Audusse and Perthame, for which uniqueness is shown 
in [12 [7D]. Thus it is not surprising that Gap is an L^D germ; property is 
checked directly from the definition (4.5) of Gap- 

Proposition 4.9. The germ Gap *s complete. 

Proof. Denote by h'i^ the inverses of /'''^ I (_oo «'•''] ' denote by h'',^ the inverses of 
/'''"Ij !,r y The domains of 6*' ■''(•, u'^'') are easy to calculate: 



ju'}, ifM'>u;,; 



Dom6''(-,u') 
Dom 6*'' (•,«'■) 



{u"-} , if u'- < <, 

^[hLiriun),^), if 

According to Remark |4.1[ we have to show that M x M is a subset of 

IJ Dom6l'(-,w') X Dom6l'^(-,u''). 

Indeed, take (u_,u+) G M x M. We have five cases: 

• in the case u_ < uj,, u+ > u^, (u_,u+) £ Dom 9^ {■ , u''^) x Dom 9"^ {■ , u^) , 
and (Mq, Uq) e Gap (notice that f'^{u''f) — 0, by the assumptions on /''''); 

• in the case U- < mJ,, u+ < u^, set s — f^{u+); then the pair (u_,it+) 
belongs to Dom 6''(-, /i'„(s)) x Dom 6'''(-, and (/i'_(s),m+) € Gap', 

• the same is true if 7i_ > uj,, u+ < and /'(m_) < f'^{u^) — s; 

• in the case u_ > m^, > m'JJ, set s — f\u-)] then the pair (u_,it_|_) 
belongs to Dom6''(-,u_) x Dom6'''(-, /i'^(s)), and (w-, /i!j^(s)) € Gap] 

• the same is true if 7i_ > Uq, m+ < u^o * ~ /'(""-) ^ /'^("+)- 

Hence, (it_,M+) e Dom0'(-,u') x Dom 6'''(-, u'') for some (u^u^) £ Gap- O 
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Notice that Propositions 4.9 and 4.51) imply that the maximal L^D extension 
(Gap)* of the germ Qap coincides with its closure Gap- The above analysis gives 

Gap - Gap U | = ^ K), > ul < <} 

EE {iu\un I f{u') = fiun, sign-(ii' - 7ii)sign+K - <) < O} . 

Now, consider the germ G[u' ,u^) ■— {{uo,u^)}- We find that {G{u' .ur))* — Gap, 
so {G(ul„u^))* is an L^D germ; by Proposition 3. 4[v), G(u'^,u^) is a definite germ. 
By Remark 3.12'ii), tj(-„i ^^r^-entropy solutions of (1.7) coincide with t^^ip-entropy 



solutions, thus they coincide with entropy solutions in the sense of Au dusse and 
Perthame [T^. The completeness of {G(u' ,u'-)) follows from Propositi on |4 .9[ The 
existence of „, )-entropy solutions can be obtained from Theorem |b.4| or from 
the adapted viscosity method of Theorem |6.3| in Section [6] (see also |70j ). 

Moreover, different notions of entropy solution can be considered for the same 
configuration of fluxes. Indeed, take any pair {A,B) e {—oo,u[] x [uJJ,oo) which 
fulfills the Rankine-Hugoniot condition; set 

■Ha,b) ■■= f{A) - r{B). 
The germ G{a,b) ■— {{A, B)} turns out to be definite, and its unique maximal 
extension {G(a,b))* is given by 

iu',U-)&iG^A,B)r 

^ f{u') = f{u-) > s^A^B), sign-(zi' - A) sign+(^.'^ - S) < 0. 
In the same way as in Proposition ■ 



(4.6) 



4.9 



we check that is complete. Clearly, 

the choice {A,B) = (Mq,Mo) discussed above is a particular case. By analogy with 
the case considered by Adimurthi, Mishra, and Veerappa Gowda in [3] (see also 
[2T] and Section 4.8 1, such pairs {A,B) will be called (A, B)— connections. 

We easily check that (4.6) accounts for all maximal L^D germs, for the fluxes 
considered in this paragraph. The monotonicity condition on can be relaxed, 
allowing for non-strictly monotone on (— cxd,^^''] and on [u^' , cxd) fluxes /'''". 

We see that the definition of entropy solutions due to Audusse and Perthame 
[T^ . as applied to the model case ( |1.7[ ), corresponds to one among infinitely many 
choices of L^-contractive semigroups of solutions. Each semigroup is most concisely 
determined by a choice of an {A, i?)— connection, i.e., by a singleton {{A, B)} from 
the set (— oo, wj,] x [u^, oo); any such singleton is a definite germ with complete dual. 

4.7. The Karlsen-Risebro- Towers entropy condition. In [7S1 [771 HSl US] , the 

definition of an entropy solution is similar to Definition 3.15 with (3.8) replaced 
by the inequality 



{\u-c\^t + q{x,u,c)^x} 



(4.7) 



\ric)~fic)\am>o, 



that should hold for aU c € M and for aU ( e 2?([0,oo) x M), ^ > 0. Uniqueness 
for (1.7) under this admissibility condition is shown for fluxes /''^ satisfying the 



so-called crossing condition, cf. Figure |2] 

3u^ £ M such that sign(z - u^){r{z) - f\z)) >0\fzeU. (4.8) 

Let us analyze the underlying germ. For the sake of simplicity, assume that 
are non-constant on non-degenerate intervals. Then the singular mappings 



are 



strictly monotone. In this case, for all u satisfying (4.7) there exist strong traces 



1 1—^ (7''''u)(t); in addition, the Rankine-Hugoniot condition holds for these left and 
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Figure 2. Left: Crossing condition (4.8) satisfied. Right: Cross- 
ing condition violated. 



right traces. As in the proof of Theorem 3.18 for > and f G I?([0, oo) x M), 

^ > 0, we take iph 
final result reads 



{h-\x\y 



as the test function in (4.7) and send h ^ 0. The 



(g'((V^)(t),c) - q^{{Yu){t),c) + \r{c) - f{c)\)at,0) > 0. (4.9) 



Setting u '^{t) :— (7 the inequalities (4.9) yield the condition 



VceM, q\{u^){t),c)~q''{{u'){t),c) + \r\c) - f{c)\>0, for a.e. i > 0. (4.10) 
A case study shows that (4.10) defines the following germ: 
(u', u'') e Qkrt if and only if f{v}) 
I either — 



and < or u' < and Vz € [u', u'"], max{/'(z), /''(z)} > s; 
[or > u'' and Vz e min{/'(z), /''(z)} < s. 



(4.11) 



This definition of Qkrt is valid also if /'''' are allowed to degenerate (become 
constant on intervals). The germ Qkrt is closed. 

Proceeding as in the proof of Theorem |3.18[ we have 



Remark 4.10. A function u is an entropy solution of (1.7 1 in the sense of |48) if 



and only if it is a solution in the sense of Definition 3.8 with Q* replaced by Qkrt- 
Thus, any function of the form (RPb-sol.) with intermediate states {u^u^) G Qkrt 



is a solution of the Riemann problem (1.7), (1.8) in the sense of [481 



Let us show that the crossing condition (4.8) is sufficient for the germ Qkrt to 
be L^D and, moreover, maximal L^D. 

Proposition 4.11. 

(i) The germ Qkrt admits no non-trivial L^D extension. 

(ii) If the fluxes satisfy the crossing condition (4.8), then the germ Qkrt is 
L^D (and therefore it is maximal L^D). 

Proof, (i) Assume that there exists {u\u^) ^ Qkrt such that Qkrt U {{u\u^)} 
is L^D. We can assume that < ii^ and that there exists z G (u^if) such that 
both /'(z) and f^{z) are strictly smaller than s /'''"(u'^'') (the other case is 
symmetric). We claim the 

3(u', u'') £ Qkrt such that < < u'' < u'' and /''''(u'''') < s. (4.12) 

Indeed, we can assume that /'(z) < /''(z) < s (the other case is symmetric). We 
set s := f^{z) and examine the three possibilities. 

- Either have a crossing point G [u\ z] such that /'■''(u^) < s. Then 



it is clear that u < < z < vT , so that u :— satisfies(4.12 1. 
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Or/'- > /'on[u',z]. Set u'^"- :=inf {C e | /''''(C) = s}. Since /''(z) = 

s, u"^ is well defined and < z < if . Since f{u'') — s > s > f\z), 
is well defined and < uK Because > /' and by the choice of m'''', we 
have u' < u\ /'■(C) > s = /'■''(u'''') VC e [u',m'']. Thus («',«'') G Gkbt 



and (4.12) holds. 



Or else, /'''' have a crossing point S [u\z\ with /'■''(■Li^) > s. Then we 
choose the closest to z crossing point with the above property, and set 
{C G [^x'-^] I /''"^(C) — As in the previous case, we check that 
are well defined, moreover, < and /''(C) > s = /''''(u'-'') for all 



Ce 



Thus, (u',m'') e feijT and (4.12) holds 



Now the pair (m',-u''), and the pair (m',m'') which satisfies (4.12)), both violate 
( |1.11[ ); thus Gkrt U {(u',u'')} is not L^D. This contradiction shows that Gkrt 
has no non-trivial L^D extension. 

(ii) Let (u^u^), (u'',u^) be two pairs satisfying the Rankine-Hugoniot condition; 
set s :— /'■'^{u''-^), s := /''''(u''''). Upon exchanging the roles of m''*" and u'''', we 
easily find that ( |1.11[ ) fails if and only if we have simultaneously s < s 
and < . These conditions are realized if and only if s < s and 



either ii < u < u"" < if 



or u'' < vJ" < < ii^ 



or u < < < u , or vJ' < u < iJ' < u , or u' 



r u'' < < u' < u'' 
< u'" < < u'. 



(4.13) 



"We check that in all these cases except for the first one and the last one, (4.11) 
never holds for both (u',u'') and (u',u''). In the two cases that were excepted. 



(4.11) can hold for both (u^u^), {u'',u^) only if the crossing condition (4.8) fails. 
This shows that if (|48l holds, then ([lIT|) holds for all («','«''), (u',m'') G Gkrt- 



We conclude that Gkrt is an L^D germ; by (i), it is maximal L^D 



□ 



Now let us look at on e case where the crossing_condition (4.8 1 fails. Consider 

uj, < u^. W e claim that 
Gkrt is not an L^D germ. Indeed, with the notation of Section [4~6| let 



the fluxes in Section 4.6 (i.e., of the type in Figure [l]) with 

!ed, with the notat 

(5 = inf {s I 3u e (Uq, Wq) such that /'(u) = /''(u) = s} ; 

then for all s e (0,(5], the pair (u^u^) := (h'-_(s), h'^{s)) belongs to Gkrt- Also 
{u\u'^) := (UgjUo) e Gkrt, but then ( |1.11[ ) is violated. Another way to prove 
that Gkrt is not L^D is to notice that the definite germ G{u' ,u^) — {{Uo,ul)} is 



contained within Gkrt, but, according to ( |4.6[ ), (u^u^) does not belong to the 
unique maximal L^D extension {G{u' of G{u' .u^)- 

Based on this fact, we can construct a non-uniqueness example for one simple 
case where the crossing condition (14. 81) fails. 



Example 4.12. Consider the fluxes f 



l.r 



Z 



\z — Uq'^'I with u[j < 



Then 



the entropy condition of 48J is not sufficient to single out a unique solution to 



the Riemann problem (1.7), (1.8) with (u_,u^) = (uq — s,Uq + s), for all s in 
(O, ^ (Mq — Mq)] ■ Indeed, the result follows by the above analysis and by Remark 
4.10 We have (u_,u+) G Gkrt; thus the function u(t,x) = U- H^^j^q} +''^+ ll{a;>o} 
is one solution. Because u^'^ S Dom0''''(-, uj,'''), there exists another solution of the 
form (RPb-sol.) with intermediate states (Uq,Uq) G Gk 



RT- 



In contrast, if we consider the fluxes as in Section 4.6 with > such that 
the crossing condition (4.8) holds, then there exists a unique crossing point u-,^ of 
/' and /'', and € (Uq,Wq). In this case, one easily checks that Gkrt coincides 
with the complete maximal L^D germ {{G{u^^u^)) ■ Another very similar example 
is given in Section [48l 
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Let US also look at the configuration of fluxes considered in Section 4.4 In this 
case, irrespective of the crossing condition (4.8 1, we always have Gkrt = Grh', 
thus Gkrt is a complete maximal L^D germ. Also in the case /' = f^, Gkrt is a 
complete maximal L^D germ, because it coincides with the Volpert-Kruzhkov germ 
GvKr (here, the crossing condition is satisfied). 

Finally, we point out that the germ Gkrt is closely related to the "vanishing 
viscosity" germ Gvv which we consider in Section [Sj more precisely, we prove there 
that Gvv — Gkrt whenever the crossing condition ( |4.8[ ) holds. 

In conclusion, the definition (4.7) of an entropy solution given in |48| can be 



equivalently reformulated in terms of the germ Gkrt (see Remark 4.10 ). Whenever 
the crossing condition (4.8) holds, the germ Gkrt is maximal L^D (sec Proposition 



4.11) and the unique entropy solution in the sense (4.7) coincides with the unique 

48] 



feflT-entropy solution (see 08] or Theorem 3.11 for the uniqueness proof). For 
general fluxes, uniqueness may fail, as in Example |4.12[ nevertheless, the crossing 



condition (4.8) is not necessary for the uniqueness of ^JxflT-entropy solutions. 



4.8. Bell-shaped fluxes. Many works have been devoted to the case of bell- 
shaped fluxes. To fix the ideas, 

let U = [0, 1]; assume f^'' : M+ are such that /''"^(O) = = /'^''(l), 

are non-decreasing on [0, u^'*"] and non-increasing on [vl'f^ 1], 

the graphs of have at most one crossing point with abscissa in (0, 1), 

and if there is a crossing point u^, then it lies between the points Mq, u^. 

(4.14) 

We refer to [211 Figure 1.1] (cf. also Figure[2]in the present paper) for a catalogue 
of different configurations. More general cases with multiple crossing points and 
oscillating fluxes are considered in . 

Here the situation is quite similar to the one described in Section [4.6| (we just have 
to take into account the fact that v}f are now points of maximum (not minimum) 
for There exist infinitely many complete maximal L^D germs. Each of these 
germs contains one and only one pair {A,B) G [wj,, 1] x [O,ito] (cf. Proposition 4.6). 
The germ is thus uniquely determined by a choice of the connection {A, B) (see 
Adimurthi, Mishra, Veerappa Gowda [3 and Biirger, Karlsen, Towers [H]); each 
maximal L^D germ is the dual {G(a.b))* of the definite germ 

GiA,B) {{A,B): A e K,l], B e [0,<], and f{A) = r{B) =: s^am)-] 
We have 

f{u') = rW) < s^A,B), sign+(7/' - A) sign- (7.'- -B)<0. 

As in Proposition 4.9 one shows that {G(^a.b))* is complete. Another way is 
to use Proposition 3.23 together with the existence result of [21]. Indeed, in view 
of Proposition 3.17 it is clear that Definition 3.15 of G(^^B)-entropy solution is 
equivalent to the definition given by Biirger, Karlsen, Towers in [21) . which requires 
that inequality ( |3.8[ ) holds 

— for aU ^ e P([0, (X)) x M), ^ > 0, such that ^(0, t) = 0; with (c', c'') arbitrary; 

— for aU ^ g ^([O, oo) x M), ^ > 0, and (c', c'') = {A, B). 

The technical genuine nonlinearity restriction on /'■'' imposed in '21] can be 
bypassecQ In Section |6j we give alternative proofs of existence of G(yi.B)-entropy 
solutions, using an adapted viscosity approach and a specific numerical scheme. 



(4.15) 



In works on conservation laws with interface coupling, such as the discontinuous flux problems, 
there are two reasons to ask for the genuine nonlinearity of the fluxes. Firstly, this assumption 
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Let us point out three specific choices of a connection {A,B). The first one is 
{A,B) = (1,0), which corresponds to the extremal choice s^^ b) = 0. Further- 
more, we have that G{i,o) = {(1)0); (0,0), (1, 1)} (this example is mentioned, e.g., 
by Panov in [TO]). The other extremal choice is to pick {A, B) corresponding to the 
maximum of S(^a,b)\ in this case, we have either A = u^, (whenever /'(wj,) < f^{u^)) 
01 B = ul (whenever /'(Mq) > /''(Mo)|3 "^^^^ germ identified by Kaasschi- 

etter [44], Ostrov [62], and Adimurthi and Veerappa Gowda [1]; the corresponding 
solutions (which are 5(A,B)-eiitropy solutions in our vocabulary) have been con- 
structed as limits of various physically relevant regularizations. In the case where 

and have no crossing point with abscissa in (0, 1), the definition of Bachmann 
and Vovelle [13] also leads to this germ, as well as the one of Karlsen, Risebro, 



Towers [48] (see Section 4.7). 



Finally, consider the case where the crossing condition (4.8) holds with S 
(0,1). Then (u^,m^) S [uq, 1] x [0,Uq]. The germ G(ux,u^) has a unique maximal 
L^D extension {G(u^.u^))* which coincides with the germ Gkrt described in Section 



4.7 As it was demonstrated by Adimurthi, Mishra, Veerappa Gowda in 3J, this 
germ also comes out of the vanishing viscosity approaches of Gimse and Risebro 
[371138] and Diehl [29l EO] [31] ■ Section [s] below is devoted to a detailed study of the 
"vanishing viscosity" germ(s). 

4.9. The case /' = and the Colombo-Goatin germs. In [25], Colombo and 



Goatin consider the conservation law of the form (1.7) with f = f = , which 



is a bell-shaped function in the sense of Section 4.8 for some F e [0, max/], the 
restriction on the value of the flux, |/(w(t, 0))| < F, is imposed. Starting from 
a regularization of the constraint and using the entropy formulation of Karlsen, 
Risebro, Towers [3S1, the authors obtain a definition of entropy solution which 
leads to well-posedness in the BV framework. 

Analyzing the definition of [26] . we easily determine the pairs (u^u^) for which 
the elementary solution ( |1.10 ) is an entropy solution in the sense of [26]. The set 



of all these pairs is in fact a maximal L^D germ; moreover, this germ is the unique 
maximal extension {Gcg) of the definite "Colombo-Goatin" germ Gcg'- 

= {{A, B): A, B are determined by ^ < S, f{A) = f{B) = F} . 

This can be seen as an (A, i3)-connection of [21] and Section |4.8[ except that we 
are in the particular "degenerate" situation with /' = /'". 

In fact, in [25] the constraint F depends on t. This case can be included in our 



consideration. We refer to Remark 3.29 and the forthcoming papers [HIIH]- 



5. The vanishing viscosity germ 
There are many possibilities for introducing a vanishing viscosity regularization 



of (1.7). Gimse and Risebro PT", '3B^ recast (1.7) as a 2 x 2 system and add the 
simplest viscosity terms to both equations. Diehl [211 [Ml HI] introduces both the 
viscosity euxx and a smoothing of the discontinuity in f; the same approach 
is pursued in ^1] . The regularization term of Kaasschietter [44] is more involved. 



simplifies the manipulation of the interface fluxes (cf. Remarks 1 2 . 2 1 [231 where these restrictions are 
bypassed). Secondly, this assumption yields strong compactness properties (see I56| and 1661 1691 ). 
In |21l . the compactness is achieved with the help of local BV estimates; thus the nonlinearity 
assumption on can be dropped. 

''We autorize to be constant on non-degenerate intervals, therefore are not uniquely 
defined. But all the germs Qf^^i (or Q(^a,u^)) corresponding to difl'erent choices of have the 
same closure. 
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making appear a capillarity pressure function; the admissibility condition he obtains 
is different from the ones of [311 [Ml EH [Ml EI] . 



We shall first study the pure vanishing viscosity regularization of (|1.7|, i.e., 
ut + f{x, u) 



eUrr 



j{x,z) 



(5.1) 



We pursue the traveling- wave approach of Gelfand [3Bj (cf. Section 4.3 1, which is, 
in our context, a "standing- wave" approach. Given /'■'', define the "standing-wave 
vanishing viscosity" germ by 

there exists a function : M M such that 

(5.2) 



lim W{£,) 



lim ly(^) = vJ" , and 



u''{t,x) = W{x/e) solves (5.1) in 2?'(E+ x M). 



Clearly, if (u ,u^) S Gyy, then the elementary stationary solution (1.10) can be 
obtained as the almost everywhere and Lj^^ limit of a sequence of solutions of the 
regularized equation (5.1). 

It turns out that Gyy can be replaced by a smaller germ which is described 
explicitly; moreover, we show that Gyy is definite. The proposition below leads 
to a definition of the maximal L^D germ Gw, which has already appeared in the 
work of Diehl [M] ■ 

Proposition 5.1. Consider the germ Gyy defined as the set of all the pairs {u\ u^) 
such that 



and 



either = u^; 
or u' < and 



or u' > ana 



either f^{z) > s for all z G (u',u''], 
or f'^{z) > s for all z G [u\ u^); 
either f^{z) < s for all z G [u'",u'), 
or f'^{z) < s for all z G {u^ , u']. 



(5.3) 



(i) The germ Gyy is L^D. 



(a) The closure Gyy of Gyy is a maximal L^D germ, and {u\u^) G Gyy 

f\u')^f-{u-)^:s, and 

' either = u^; 
or u' < and there exists u° G 

such that f^{z) > s for all z G [u', u°] 
and f^{z) > s for all z G [u°,u^]] 
or > u'^ and there exists u° G [w' , u'] 

such that f^{z) < s for all z G , u°\ 
and f\z) < s for all z G 



(5.4) 



(Hi) We have Gyy C Gyy C Gy 



'vv 



(iv) Both Gyy and Gyy are definite germs with the same unique maximal L^D 
extension Gyy- The Gyy, Gyy, and Gyy entropy solutions coincide. 

Definition 5.2. In the sequel, we denote by Gvv the maximal L^D germ Gy 



vv 



given by (5.4). We call it the "vanishing viscosity" germ associated with the pair 
of functions /''''. 
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Whenever the dependence on f'-'^ should be indicated exphcitly, we denote the 
corresponding vanishing viscosity germ by Qyy ■ In [9] , we will study more carefully 
the dependence of Gyy on /''^. 

Remark 5.3. 

(i) Definition 5.2 allows for general continuous fluxes Yet it is clear that the 
so defined germ Gvv need not be complete. For instance, if the ranges of do 
not intersect, Gvv is empty and no Gvv entropy solution exists. 

(ii) Basing themselves on different viscosity approaches, Gimse and Risebro, and 
then Diehl, have implicitly introduced the set of all admissible elementary solutions 
of (1.7). These are the elementary solutions used to construct a Riemann solver 
based on the "minimal jump condition" |371I38| or on the 'T -condition" [291 1301 15T] . 
In fact, both approaches lead to the same germ which coincides with Gvv j the closure 
operation accounts for the difference between the minimal jump condition and the 
T-condition. For certain configurations of the fluxes /'■'', a Riemann solver was 
defined for all data in [371 |3H1 HH EHl [SI 132] • o-H these cases we conclude that 
Gvv is complete. 

(Hi) In [32 , Diehl reformulated the T-condition from his previous works [2911301 [5T| 
in the following form: the pair (u', u"^) satisfies the T-condition if 



f\u^) — f^{u^) and there exists u° G ch{u^,u^) such that 
(u'' - u°) ( fiz) - > 0, Vz e c/l(^^^ u°) 



(5.5) 



where for a,b G M, ch{a,b) denotes the interval [min{a, 6}, max{a, 6}] . Clearly, 



(5.4) coincides with (5.5). The descriptions (5.4) and (5.5) are reminiscent of the 
Olemik admissibility condition (for convex fluxes f'"^) and of the "chord condition" 



(see [33 1301 [31] as is the case with the chord condition, (5.4) and (5.5) are derived 
from the travelling-waves approach of |36j . 
Proof of Proposition |5.1[ 



(i) The proof is entirely similar to the one of Proposition 4.11 u). Because / 



are continuous, each of the six cases in (4.13) turns out to be incompatible with 
the assumption that both {u\u^) and {u^,u^) fulfill (5.3). 



fl.r 



Notice that a less tedious proof can be obtained by the passage to the limit as 



e — >■ in the Kato inequality which holds for equation (5.1) (indeed, recall that 
( 1.11 ) follows from the Kato inequality ( |1.5[ )). Such a proof is given in ^TU[, for the 
case of Lipschitz continuous /''''. The general case follows by approximation. 



(ii) Recall that we denote by Gvv the set of all pairs satisfying (5.4). In a first 
step, we show that Gvv C Gyy- Since Gyy is an L^D germ by (i) and Proposition 



4.4 i), this implies that Gvv is an L^D germ. 

Then in a second step, we show that Gvv contains the dual (Gyy) 



of GZ 



According to Propositions 



4.4 



and 



3.4 



vv 



this yields the reverse inclusion Gvv ^ Gy 



and hence the maximality of the L^D germ Gvv — Gvv 

We will make repeated use of the continuity of /''^, without mentioning it. 



vv 



Step 1. Let (u ,u^) e Gvv- If w ~ u'', then [u ,u^) G Gvv Gvv ^^'^ remaining 
cases are symmetric; let us treat the one with < . 



Take s and u" as introduced in (5.4). Consider the function / (•) — s on the 



interval [u°,u'^). If it has a zero point, set z' :— minjz e [u°,u'^) | /'(z) — s} 
and := min {z £ [z\u^] | f^{z) = s} (since f^{u^) = s, z^ is well defined). By 
construction and by (5.4), we have 

— < z', f^u^) = /'(z') = s, and /' > s on [m',z']; 

— z' < z*", /'(z') = f^{z^') = s, and f^ > s on the interval (z', z"^); 
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— z'' < M^ r{z'') = /'■(«'■) = s, and f > s on [z'', u'']. 

This means that (m',z') (resp., {z'^,u^)) is a left contact shock (resp., a right 
contact shock), and {z\ z^) e Qyy 

Now consider the situation in which /'(•) — s has no zero point on the interval 
[u°,u'''). Then /' > s on [u^u^]. In the case f{u'^) = s, {u'-,u^) is a left contact 
shock and (u^u^) € Gyv Otherwise, we define z' := max {z e [w', w'') | f\z) = s}. 
In this case, (w',z') is a left contact shock. In addition, (z',u'') e CJyy because 
z' < u"^, f\z^) = /''("'') = s, and /' > s on (z 



In all the cases, by definition of the closure, we conclude that (u', u^) S Qyy- 
Step 2. It sufhces to show that if /'(u') — f^{u^) but (u^u^) ^ Gvv, then there 
exists (c^c^) e tyyy such that 

q\u\c')<q'-iu\c^). (5.6) 
Set s /''''(u''''). As before, it suffices to consider the case u' < u^. Define 
z' sup {z e [u', u'"] I /' > s on [u', z]} , z*" := inf {z e [u', u''] | /'' > s on [z, u'']} . 



If z > z*", then (5.4) holds with u° — z , so that (u , u'') lies in Gyy. Thus, z < z^. 
Now, there are three cases to be investigated: 

(a) /' and have a crossing point z" in the interval {z\ z^) such that /''''(z°) < s; 

(b) /' and have a crossing point z° in the interval (z', z'') such that /''''(z°) > s; 

(c) either /' < on the interval (z', z''), or Z*" < /' on the interval (z', z^). 
J 



In case (a), setting (c , c*") :— {z°,z°) we obtain ( |5.6[ ), because m < c , c'' < u'' 

o 

VV 



and /'''■(«'''■) > /'■'■(z°). Notice that we do have (z°,z°) € 



In case (b), by definition of z''^, there exists s < s such that s belongs to 
/' ([z',z°]) n /'■([z°,z'']). In this case, set c' := max{z e [z',z°] | /'(z) = s} and 
:= min {z £ [z°, z''] | /''(z) = s}. We then have ( |5.6| for the very same reasons 
as in case (a). In addition, (c',c'') G ^Jyy, because < c"" and /' > s on (c',z°], 
/•^ > s on [z°,c''). 

In case (c), the two situations are similar. Consider, e.g., the case /' < on 
(z', z^). Choose for the point of [z\ z^] where attains its minimum value over 
[z', z"^] and which is the closest one to z'. By definition of z*", s := f^{d') is smaller 
than s. Because f^{z'') > f^{z^) — s > s, we have > z'; in turn, this yields 
f''{c^) < f^id") — s. Since f^{z^) — s > s, there exists c' in the interval (z^z^) 
such that f{c^)^s. The pair (c', d") fulfills ( |5.4| . In addition, by the definition of 
d we have Z'' > s on [z\ z^] D [c', c""]; thus ((tTc^) G ^yy. 

In all cases, we have constructed (c',c'') S Qyy with property (5.6). 

The contradiction shows that {Qyy) C Qyy. Thus (ii) follows. 

(iii) We first show the inclusion Qyy C Qyy Let {u\u^) satisfy (5.3). In the 
case li' — u'' , the standing- wave profile W can be chosen constant on ]R. The four 
other cases are symmetric. For instance, in the case u' < and /'(z) > s for all 
z € (u',u''], is a continuous function that is constant (equal to u*") on [0, oo). 
On the interval (— oo, 0], is constructed as a solution of the autonomous ODE 

W' = f{W)- f{v}), W{Q)=vr. 

Indeed, because f{W) - /'(«') = f{W) - s > for taking values in (u', u''], W 
is non-decreasing. If the graph of W crosses the line W — v} extend the solution 
by the constant value v} on the left from the crossing point. Otherwise, W is defined 
on the whole interval (— oo,0], and there exists d := limj_j._oo W^(C) £ [u\u^]. In 
this case, /'(d) — s = 0, which yields d = uK Therefore (5.2) holds. 

Next, we show that Qyy C Qyy Let {u\u^') e Qyy The case = is trivial; 
by a symmetry argument, we can assume u' < u'' . The standing- wave profile W in 
(5.2) is a continuous on M function which is monotone on M_ and on Moreover, 
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replacing W by max {u', min {VK, it' }} , we still have a standing- wave profile with 
the same properties. Then u° := W{0) lies within [u', u^]. The monotonicity of W 
guarantees that — s > on [u', and /'' — s > on [u°, u^]. Thus w' ) fits 
the definition Q. By (ii), (u',it'') e G^. 



(iv) The claim follows readily from (ii) , (iii) , and Proposition |4.4| (iii) , (iv) . □ 

Now we compare Qvv with the known germs studied in Section |4] and in |37| 
l38j l29l [30l [31] l32] . Let us emphasize that in many cases studied in these papers, 
Gvv turns out to be complete. The completeness of Gvv will be further studied in 
Theorem |6.1| (see Section [6]) . 



Remark 5.4. 

(i) With the help of the characterizations (4.11), (5.3), (5.4) oJQkrt, Qyv' ^'^'^ 
Gvv J respectively, we easily see that Qvv Gkrt = Qvv whenever the crossing 
condition (4.8) holds. In this case, the Gvv md Gkrt entropy solutions coincide. 



(ii) Consider the case of bell- shaped fluxes (4.14) treated in Section 4. 



— First assume < Ug and that there exists a crossing point G 
between /' and f^ (in this case, the crossing condition (4.8) holds). Then Gvv 
contains the definite germ G(u^,u^)- Thus the G{u,^,u^), Gvv, o,nd Gkrt entropy 
solutions coincide. In addition, since {G{u^^u^)Y is complete, Gvv is complete. 

— //, on the other hand, a crossing point exists but u^^ < Uy, < (this is 
the "bad" crossing case), then the Gvv -entropy solutions coincide with the entropy 
solutions for the germ G(a,b) with {A, B) corresponding to S(a,s) — Smax, where 

w:-niax{seM+ | 3(A,B)e K,l] x [0, <] s.t. f\A) = f^ [B) = s] . (5.7) 

Also in the case where the fluxes do not cross within (0,1), Gvv = G(a.b) with 
^{A,B) = Smax- both cases, Gvv 'is complete because {G(a.b)) is complete. Let 
us stress that in the non-crossing case, Gkrt coincides with Gvv , whereas in the 
case of a "bad" crossing, Gkrt \ Gvv 7^ 0- 

In conclusion, following Diehl [5^11301131] and Biirger, Karlsen, Towers let 
us mention that the germ Gvv ^^^^ ^-rises from the more general regularization- 
viscosity approach 

ut + f {u)^ = en.,, in 2?' (M+ X M) , (5.8) 
where is a continuous (in x) approximation to f. 

Proposition 5.5. Suppose /'''' are locally Lipschitz continuous. Let (f'^(-, •))(5>o ^6 
a family of continuous functions approximating f in (1.7) in the following sense: 

- f{x, z) is of the form F{x/S, z) with some fixed function F; 

- for all y e K, there exists ay e [0, 1] such that F{y, •) = ciyf'-{-) + (1 — 

- for y < -1, ay ^ I (i.e., F{y, •) = f{-)); similarly, for y > I, F{y, •) = 

n-)- 

Let (u^u^) G Gvv Denote by L the Lipschitz constant of f'"^' on the interval 
between u'' and . 

Then, whenever 5 / e < 2/L, there exists a stationary solution u{t,x) — u^'^{x) of 
(5.8) satisfying lima;_>._co u^'^{x) — w', lima;_>.co u'^'^ix) — . In addition, as e J, 0, 
with < (5 < j^e, u^'^ converges a.e. and in L\^^{U.) to the elementary stationary 



so. 



Zwizon 'u']l{^<o} + ^^''ll{a;>o} of ( [TJ ). 

o 

VV 



Proof. Let {u^u^) G Gyv "^^^ ~ "^^ trivial since constants solve (5.8). 



The other cases are symmetric, so let us only consider the one with < and 
/' > s on the interval (u' , u^] . 
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Setting y — u'^'^{x) — U{x/S), we are reduced to the problem 

e U'iy) = S {F{y, U{y)) - s), U{~^) = u\ U{^) = . (5.9) 

In view of the assmiiptions on F and also the ratio 5/e, the (constant) function 
m+ : y G (—00, 1] i-)- u*" is a subsolution of the ODE U' ~ 6/e{F{y, U) — s); and 
the function m_ : y £ (—00, 1] i— > min{w',M'" + L5/e [y — 1)} is a supersolution of 
the same ODE. By the classical ODE existence and comparison result, there exists 
a solution of the ODE on (—00, 1] with U{1) = u"", and we have m_ < U < m+. 
In particular, U{-) takes values in {u\u^] and satisfies the autonomous equation 
U' = S/e{f'-{U) — s) on (— 00,— 1). Since by assumption, /' > s on {u^u"^], U{-) 
admits a limit as ?/ — )• —00 which necessarily equals uK Finally, extending U by the 
constant value on [l,oo), we obtain a solution to problem ( |5.9[ ). 

For the proof of convergence of the solutions u^'^{x) — U {x/5), we replace m+(-) 
on (— 00,— 1] by the solution V{-) of V = 5/e{f{V) — s) with the initial datum 
V{—1) — vT] with the same arguments as above, such a solution exists, it tends 
to as 2/ — > — cx), and it bounds t/(-) from above. By assumptions, we have 
6 < 2/Le as e — 0. For all a; < 0, for all sufficiently small e, the quantity 
\u^'^{x) — Ti'l is upper bounded by \V{x/6) — w'| which converges to zero. For 
a; > 0, we simply have u'^'^{x) — whenever S < x. We conclude the pointwise 
convergence of u^-^ to the profile (|1.10|); since all these functions take values within 



the Lj^^ convergence follows. □ 



6. Some existence and convergence results 

In this section, for a fixed germ Q, we study the existence of C/-entropy solutions 
and convergence of various vanishing viscosity and numerical approximations. 

The first result concerns the vanishing viscosity germ Gvv, which was described 
in Section [sj The main assumption for the existence result is the availability of a 



uniform L°°(IR+ x M) estimate for solutions of (5.1 1 with a given initial datum 
Mo G L°°(IR). We also impose some additional nonlinearity and Lipschitz continuity 
assumptions on the fluxes /''''. For the sake of simplicity, we restrict our attention 



to fluxes satisfying /''""(O) — — f'''^{l) (cf. the bell-shaped fluxes (4.14)) and an 
initial condition uq with values in [0, 1]; thus solutions are automatically bounded, 
because they take values in [0, 1]. 

Recall that in the setting ( |4.14 ) of bell- shaped fluxes with < < u^, the 
germ Qyv coincides with {G{u^,u^))* the other maximal L^D germs {G{a,b))* 
are corresponding to other choices of the connection {A,B). For an arbitrary but 
fixed connection {A,B), using a specially fabricated artificial viscosity we prove 
convergence of viscosity approximations to the C/(^_B)-entropy solution. 

Then we look at a general complete L^D germ G for a pair of locally Lipschitz 
continuous functions /''''. We construct solutions by showing convergence of a 
suitably adapted monotone three-point finite volume scheme. For these results, 
the uniform BVioc estimates away from the interface {x = 0} are of importance 
(see [inilSI])- As previously, we need a uniform L°° bound on the approximations. 



which actually comes for free from the completeness of G (cf. Proposition 3.201. 

Finally, we prove convergence of uniformly L°° bounded (thus weakly compact) 
sequences of approximate ^-entropy solutions, without utilizing BV type a priori 
estimates and assumptions that ensure strong compactness. For this purpose, we 
have to assume that C/ is a maximal L^D germ for which existence is already known 



(cf. Section 1.4). Then, assuming that the approximation procedure is compatible 
with the germ G (i.e., that the elementary stationary solutions selected by the 
germ G are obtained as limits of the approximation) and that an appropriate Kato 
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inequality holds for the approximate solutions, we deduce the convergence of the 
approximating procedure. 

6.1. The standard vanishing viscosity approach. To be specific, let us work 
with the vanishing viscosity germ Gvv in Definition |5.2[ and note that according 
to Proposition 5.f iv), we can replace Gvv by Gyy or by Gyy 

Theorem 6.1. Suppose : [0, 1] — > M. are Lipschitz continuous functions that 
are not affine on any interval I C [0, f]. Moreover, assume /''^(O) = = /''''(f). 
Then for all Ug which is measurable and takes values in [0, f], there exists a unique 



Gvv -entropy solution to problem (f.7),(f.2). In particular, the germGvv is definite 
and complete. In addition, the Gw-en-ti^opy solutions are the vanishing viscosity 
limits in the following sense. For each e > 0, there exists a weak solution € 



r2 



+ of the viscous problem (5.f I with any measurable initial data 

u'^lt^o — Uq taking values in [0, f]. // we moreover assume Uq — > Ug a.e. on M, then 
u"^ converge a.e. on x M <o the unique Gvv-entropy solution of ( |f.7[ ),(f.2). 



A multi-dimensional analogue of Theorem 6.f holds, see \TU[ 1^ 



Remark 6.2. 

(i) The assumptions /''' (O) — = /''' (f) o,n-d Uq{x) G [0, f] are put forward to 
ensure that the family (M^)£>g is bounded in L°° . Indeed, Theorem 6.f remains true 
if the L°° bound is provided by some other means. This boundedness assumption 
is important. For example, in situations where the graphs of /' and f^ do not 
intersect, it is then clear that (1.7) has no solution and consequently ||u^||oo oo 
as e ^ 0. 

(a) On the contrary, the Lipschitz continuity assumption on /'''" is only needed to 



establish the existence of a solution semigroup for viscous problem ( 5.f ) that satisfies 
the Kato inequality. This assumption can be bypassed, thanks to a regularization 
of /'>''. We then establish the existence of a (possibly non-unique) solution 



to (5.f ),(f.2) with merely continuous functions /'■'', and the convergence of the 
sequence {u'^)e>o to the unique Gvv -ent ropy solution of problem ( |f .?[ ) ,( f .2 ). 



(Hi) As pointed out in Corollary 



6.6 



the assumption of non-degeneracy of /'■'' 
imposed in Theorem 6.f can be dropped; but the proof becomes much more indirect, 
involving a beforehand justification of the existence of Gw^i^tfopy solutions (for 
this, we use a numerical scheme, see Theorem 6.4). Moreover, in Corollary |6.6| 
we do not prove that the germ Gvv is complete, but we take completeness as a 
hypothesis. 



Proof of Theorem 6.f Extend /'''' by zero outside [0, f]. Let T denote a generic 
positive number. The following a priori estimate is easily obtained: 




(6.1) 



r{z) dz 



where u*^ G (^(M^; iJ^(M)) is a solution of (5.f ). Thanks to this estimate, it can be 
shown by the classical Galerkin technique (see, e.g., Lions [5S]) that for any e > 
and Ug G L^(]R) there exists a solution G C(M+; 7J^(K)) to the Cauchy problem 
for the parabolic equation ( |5.f [ ). What's more, the Kato inequality holds for |5.1|). 
Indeed, assuming (ug — Ug)+ belongs to G i^(IR), let G Lf^^{R+; H^^^ 
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be two weak solutions of (5.1) corresponding to the initial data Uq,Uq S L°°(E), 
respectively. For a < 0, let us introduce a Lipschitz continuous approximation Ha 
of the sign^(-) function: 



Ha{z) := 

and use the test function (m*^ 



1 



z e 



(6.2) 



) , in the weak formulation 
of (|5.1|. Notice that, due to the Lipschitz assumption on 



<ii(/')'ik-iieii 



{0<u=-«=<a} 



iu^-u^)J — ^0 asa^O, 



thanks to (6.1); and the same holds with replaced with M-|_,/''. Hence as 



Q — > 0, we deduce the Kato inequality: for all ^ S I?([0,oo) x M), ^ > 0, 



+ / K - ■)>e / [{u^ - uT] , i.; 



(6.3) 



where we have written for the semi-Kruzhkov entropy flux: 

q+(x,z, k) := sign+(z - k){{f\z) - /(fc))ll{,<o} + {f {z) - T (fc))ll{,>o}) • 

In this inequality, using, e.g., the technique of Maliki and Toure [58] we can let ^ 
converge to the characteristic function of [0, T) x M in such a way that the terms 
with vanish as e — > 0. Then we get the contraction and comparison inequality 



for a.e. t > 0, 



(6.4) 



Inequality ( 6.4 ) ensures the uniqueness of a weak solution in L^^^ (IR+ ; Hl^^ (M) ) of 



the Cauchy problem for the parabolic equation (5.1). It also yields the comparison 



principle. Keeping in mind that, under our assumptions, the constants and 1 are 



evident solutions of (5.1), we derive the maximum principle: for data Mq taking 



values in [0, 1], < u"^ < 1 holds a.e. on 



Now we are in a position to extend the existence result for (5.1) to a general 



initial function uq taking values in [0, 1]. Indeed, following Ammar and Wittbold [S] 
we can take the bi- monotone approximation of uq € L°° (M) by bounded compactly 
supported functions, namely, 

(uo)m,n := min{(uo)+,n} 11{|2:|<„} - min{(uo)",m} ll{|j;|<™}. 

We deduce the existence of a weak solution to Moreover, thanks to the mono- 



tone convergence theorem, the obtained solutions still fulfill (6.3) and (6.4). Finally, 
the solutions belong to L^Qp(]R+; iJi^j,(]R)), thus the uniqueness result applies. 

Now, let us justify the convergence of to a tjyy-entropy solution. Take a 
sequence of data Uq as in the statement of the theorem. By the non-degeneracy 
assumption on /'''' and strong precompactness results of [521 HI], applied separately 
in the domains {x > 0} and {x < 0}, we deduce that there exists a (not labelled) 



sequence e \, such that converge to some limit u in Lj^^ 



By the 



dominated convergence theorem, we also have the L^^^ convergence of Uq to uq 
Passing to the limit in the distributional formulation of the Cauchy problem for 
(5.1), we find that m is a weak solution of (1.7), (1.2). Following [SD], we also 



deduce the Kruzhkov entropy inequalities in the domains {a; > 0} and {x < 0}. In 
particular, the existence of strong traces 7''''u on {x = 0} follows. 
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Finally, notice that for any pair u, u obtained as the vanishing viscosity limit with 
the same extracted subsequence e — 0, we can pass to the limit in the Kato in- 



in Li 



equality (6.3), where the right-hand side vanishes due to a uniform bound on £|u|P 
X M). To conclude the proof, it remains to notice that, by the definition 

c'^llr 



(5.2) of Gvvi elementary solutions u = c'llj^-^o} +c.'~^{x>a}i [c ,c^) € Gyy, are 
obtained as vanishing viscosity limits. Passing to the limit in the corresponding 



and zero remainder term TZg. According to Proposition 



3.15 



u is a CJyy-entropy solution of (1.7), (1.2). By Proposition 



Kato inequalities written for £ > 0, we infer the entropy inequalities (3.8) with 
Definition 



3.17 



5.1 



and 
tiv), 



{Gvv)* coincides with the maximal L^D germ GvV: thus we conclude that u is the 
unique Qyv-entYopy solution of il.71, (11. 21). The uniqueness of the accumulation 



point ensures that all sequences converge to the same limit u. 

It remains to notice that we have obtained a solution for every measurable initial 
function uq with values in U — [0, 1]. In particular, for all Riemann initial data uq 



in (1.8), there exists a ^yy-entropy solution. According to Remark 3.13 the germ 
Gvv is complete. □ 



6.2. The vanishing viscosity approach adapted to {A, i?)-connections. In 

this section we work with bell-shaped fluxes as defined in (4.14). Recall that any 
complete maximal L^D germ is of the form {G(^a,b)) j where G(a,b) 
a definite germ and the connection {A, B) is a pair satisfying 



:={(Ai3)}is 



{A,B) e K,l] X [0,<] and f{A) = r{B) =: s 



(A,B)- 



(6.5) 



The explicit description of (G(^ s)) is given by formula (4.15) 



Our goal is to construct B)-entropy solutions by the vanishing viscosity 
method; clearly, the choice of viscosity must be adapted to the connection (A, B). 
For the Buckley-Leverett equation with a flux that is discontinuous at a: = 0, 
Kaasschietter [44] gives a physically motivated viscosity of the form 

e{X{x,u)prXu)x)^, 

where Pc(') is the capillary pressure and \{x,-) = A'(-)ll{2,<o} + ^^{')^{x>o} is 
the mobility function, discontinuous at {x = 0}. This choice corresponds to the 



particular connection {A,B) such that S(^A g^ = Smax) cf. i5.7h 



Our approach is more academic. We fix an arbitrary connection (A, B) and 
construct an artificial "adapted viscosity" of the form 

e{a(x,u))^^ 

such that the stationary solution c{x) := Allj^^Qj + Bll^^^^Qj of the limit equation 



(|1.7P is also a solution of the viscous problem 

Ut + f{x,u)^ = ea(x, 



fix,z) 



r{z), x>o. 



a{x, z) 




(6.6) 



In the sequel, by a weak solution of (6.6) we mean a function u e L°°(M+ x M) 
satisfying (|6.6| in the sense of distributions and such that w{-) :— a{- ^u{-)) belongs 



For (6.6) to be parabolic, the functions a'-*" should be strictly 
increasing on [0, 1]; and it is convenient to ask that inf^gjo^ij (a'''') [z) > 0. Because 
/''""(O) = = /'''^(l) by assumption, it is convenient to require the constants u = 



and u = 1 to be solutions of (6.6) for all e > 
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For example, any continuous functions a ''' such that 



'a'(O) = = a''(0), 



are strictly monotone on [0, 1], and < a' (A) 



a\l) = l = a-{l), 



(6.7) 



with K e (0,1), do satisfy the above requirements, except with [A^B) — (1,0), 
in which case we take k — Q. For example, we can interpolate the values (6.7) 



to produce piecewise afhne and strictly increasing examples of such functions a^'^ 
defined on [0,1]. Notice that, by choosing k, appropriately, for the connection 
(m^,u^) we obtain the standard viscosity euxx studied in the previous section; 
under assumption (4.14), we actually have (Q{u^,u^))* — Gvv- 



Now we will prove that the ^(A._B)"6ntropy solutions of (1.71 can be obtained 
as the limit of the vanishing adapted viscosity approximations defined by ( |6.6| . 
More precisely, we have the following analogue of Theorem |6.1[ which represents a 
reinterpretation and improvement of the results of Biirger, Karlsen, Towers [2T\. 

Theorem 6.3. Assume that /'''' are Lipschitz continuous functions of the form 



(4.14). In addition, assume that /' and f^ are not affine on any interval I C [0, 1]. 
Let (A, B) be a given connection of the form (6.5). Then, for each measurable initial 
function uq taking values in [0,1], there exists a unique Q(^A.B)-£nt^opy solution of 
(1.7), (1.2). Moreover, consider a'' satisfying (6.7| and set 



a{x, ■) = a'(-)ll{^<o} + a'~(-)ll{a;>o}; 

assume in addition that {a^'^)' > const > a.e. on [0, 1]. Let (uq)£>o be a family of 
measurable functions taking values in [0, 1] such that Uq — > Mo a.e. on M. For each 
e > 0, there exists a weak solution of the adapted viscosity regularized problem 
(6.6) with initial data m^L „ — Uq- The sequence {u^)e>o converges a.e. on M_|_ x M 



to the unique G{A.B)-^''^^'''opy solution of (1.7), (1.2) 



Proof. We will focus on the existence of solutions to the adapted viscous problems 
and the convergence of these solutions as the viscosity parameter tends to zero. 
The proof is very similar to the proof of Theorem |6.1| We will argue in terms of 
the unknown w{t,x) := a{x,u(t,x)). Then 

u{t,x) ^b{x,w{t,x)), b{x,-)^{a'^) ^ (•)ll{^<o} + (a'')"^ (•)ll{2:>o}, 



and the Cauchy problem for (6.6) is equivalent to 

b{x, w)t + f{x, b{x, w))x = ewxx, w{0, x) = Wg := a{x, ul{x)). 



(6.8) 



Existence for (6.8) is obtained by proving convergence of Galerkin approximations 
along the lines of Alt and Luckhaus [1]. Indeed, set B{x,r) := b{x, z) dz. In the 



same way as for (6.1), assuming that Bo Wq{-) := B{-,Wq{-)) belongs to L^(M) and 
using the weak chain rule we obtain the a priori estimate 



\w. 



,e|2 



< 



B{;w'{T,-))+e 



Bi;u%0,-)) 



/ / f\z)a\z)dz] + i r{z)a^{z)dz 

J-oo \Jo J Jq Jo \Jo J 
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This estimate and the compactness technique of [11 ensure the convergence of the 
Galerkin approximations to a weak solution w'^ of (6.8 1 with the properties 



Bow" e L°^{R+]L\R)), w" e Ll^{R+;H\R)), 
[b o w% e LlM+\H-\R)), {b o w'){0, ■)^bo w'^. 
As in the proof of Theorem 



6.1 



the assumption S o g L^(M) can be dropped, 
thanks to the bi-monotone approximation of Wq by compactly supported bounded 
functions 

{wo)rn,n min { (wg) + , n} ll{|a;|<„} -min{(u;g)",m} 11{|:e|<„i} 

(see [5]) and thanks to the comparison principle that we now justify. Indeed, let 
w'^ be a solution of the same equation corresponding to the initial datum Wq. Take 
f £ r'([0,oo) X M), ^ > 0, and Ha given by (6.2 1. Utihzing the test function 
Ha(w" — w'^)^, using the doubling-of- variables technique in time only, as in Otto 
|64j . and taking into account that /'■'' and 6'"'' are Lipschitz continuous (cf. the 



proof of Theorem 6.1), we deduce the Kato inequality 

((& o ui^ — 6 o ■w")'^^t + b o w", b o w") 



+ / (6owg-&o{tg)+e(0,-) >e 



here q+ has the same meaning as in (6.3). With the technique of FSHI, we can let ^ 
converge to the characteristic function of [0, T) xR and derive the contraction and 
comparison inequality 



for a.e. t > 0, (u" - 11")+ (t) < (uf^ - u^) 



By the definition of {A, B) and assumptions ( |4.14 ), (6.7), the constants and 1 are 
evident solutions of (|5.1|). Hence, the following maximum principle holds: for data 

< 1 a.e. on R+ x R. Also the function 
As in the proof of 



taking values in [0, 1], we have < 



;<o} + -B ll{2;>o} is an apparent solution of (5.11 
we use uniform estimates on u' 



u{x) := A\ 

Theorem |6.1[ we use uniform estimates on and compactness arguments to pass 
to the limit in the Kato inequality written for and for = u. What we obtain 

(this is exactly the 
we see that m is a 



is the entropy inequality (3.8) with (c , d") 



{A, B) and Tig 

entropy inequality of Biirger, Karlsen, Towers [21] )■ In addition, 



weak solution of (1.7) 
{a; < 0} and {a; > 
using Proposition 



3.17 



(|1.2|) and it is a Kruzhkov entropy solution in the domains 

, or 
□ 



0}. We conclude either using the uniqueness result of 
and the fact that Q(a,b) is definite. 



6.3. Existence for complete germs through the discretization approach. 

We now establish the existence of a C/-entropy solution for any complete maximal 
L^D germ Q. For the sake of simplicity, let us take [/ = M. 

Theorem 6.4. LetQ he a complete maximal L^D germ. Assume that the junctions 
are locally Lipschitz continuous on R. Then for any initial function uq G L°°(]R) 



there exists a unique Q-entropy solution of problem (1.7|,(1.2| 



Proof. The proof is a combination of the well-known finite volume method (see 
Eymard, Gallouet, Herbin [331) with a careful treatment of the interface {x = 0} 
using the Godunov scheme; the proof of compactness of the family of discrete 
solutions is based upon the BV\oc estimate of Burger, Karlsen, Towers [THl HI] and 
an L°° bound similar to the one of Proposition 3.20 



Recall that by the definition of a complete germ, there exists a Riemann solver 
TZS'^ at {x = 0}, fully determined by CJ, defined for all Riemann data (1.8). 



We 



42 B. ANDREIANOV, K. H. KARLSEN, AND N. H. RISEBRO 

set m := essinfRUo (resp., M :— esssupgUo) and denote by c'"'^ (resp., by C'''') 
the one-sided traces "/'''^TZS^m (resp., j''''^TZS^ M) of the solution of the Riemann 
problem at {x = 0} with the constant datum = m (resp., with the constant 
datum = M). Assign 

:=min{c''^m}, S'^'' max {C'^^ M} . (6.9) 

We consider an explicit finite volume scheme based on the uniform spatial mesh 

|J.^^(iAa;, (j + l)Ax), Ax > 0, 

and a time step At satisfying the CFL condition 

A:=^<l/2i, L:=max| max \{f)'{z)\, max |(D'(z)| 

A.T \ze[VM^] ' ' z<£[b^,B^] 

A lower bound for L can be obtained from Proposition |3.20| 

For z ^ 0, we can utilize any consistent monotone flux (?i+i/2(-, •) (see e.g. [55]): 
but at the interface z = we take 5i+i/2(-, •) to be the Godunov flux based on the 
exact 5-Riemann solver. Given the numerical flux .gi+i/2(-, the difference scheme 
is defined as 

< = ur' - A - «"\ <ri')) , Vn e N, Vz e Z. 

where the iteration is initialized by 

M° := -r— / Mo, Vl G Z. 

JiAx 

The monotonicity of the numerical flux .gi+i/2(', •) and the CFL condition ensure 
that the scheme can be written as 

< = «ri\ , Vn G N, Vz G z, 

for some functions Hi{-,-,-) that are monotone in each of the three arguments. We 
identify the discrete solution (uf)neN,iez with the piecewise constant function 

° ^ * ({n-l)At,nAt)x{iAx,{i+l)Ax))' 

where /i is a collective symbol for the discretization parameters Ax, At. 
By construction (exact C/-Riemann solver at the interface {x — 0}), 



(6.10) 



the difference scheme preserves all stationary solutions of (1.7) 
of the form c' ll{2;<o} + c'' 11{:e>o} with (c',c'') G Q. 
First, let us establish the following uniform L°° bound on S'^Uq: 
b{x) -.^b' ll{:r<o}+6'' ll{x>o} 

<{S''uo)it,x) (6.11) 
< B' ll{,<o} + B'- ll{,>o} Bix), 



for a.e. (t, x) G (0, oo) x M, where b''^,B'''^ are defined in (6.9 1. To prove (6.111, 
we take the Riemann initial datum b{-) and look at the corresponding discrete 
solution {S'^b){At, ■) after the first time step. Denote by 6f the values taken by 
{S^b){At, ■). Notice that, by the choice of fe''*", the exact solution TZS^b of the 
Riemann problem has the one-sided traces c'"'" at {x = 0}. Therefore, the Godunov 
flux gi/2{bQ,bi) = gi/2{b'' jb"^) at the interface at time level n = 1 takes the value 
/'■'"(c'''') (the two values being equal). From the definition of the scheme, it is clear 
that bj = b^ for alH 7^ 0, 1. Moreover, from the definition of the scheme we have 

bi = b'-\ (/'(c') - f\b')) , b\ = b^- \{rm - nc^)) ■ 
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Because 5' < c' by construction, and because the state 6' at x < can be joined 
to the state c' at a; = 0— with waves of negative speed for the flux it foUows 
that /'(c') < /'(6')- Similarly, we have f^{c^) > f^ib"^)- Combining the above 
information, we deduce that for all i G Z, bj > . Thus we have 

{SH){At,-)>b{-). 

Hence by induction, using the monotonicity of the scheme, we deduce 

{S''b){t,x) > b{x) for a.e. {t,x) G (0,oo) x M. 

Because uq > b a.e. on M and again in view of the monotonicity of the scheme, 

{S^uo){t, x) > {S''b){t, x) > b{x) for a.e. {t, x) e (0, oo) x E, 



which establishes the desired lower bound in (6.11 1. The proof of the upper bound 
is entirely similar. 

Now we assume that uq is compactly supported and belongs to BV{M.). Then 

for afl r > and / > 0, we apply the uniform Bv(^{0,T) x (M \ (-Z,0)) estimate 

of [TH HI] to the family {S'^uo)h>o- To establish this estimate, we first combine 
the discrete contraction property, obtained from the Crandall-Tartar lemma 
[28], with the Lipschitz continuity in time t of the solver S'^. The result is 
an estimate of the time variation of S'^uq on (0, T) x K in terms of the variation 
Varzio of the initial data uq on K. Then, using the mean value theorem, we pick 

£ (0, 1) such that the variation of {S'^uo){-,±l'^) does not exceed Const (Varug) /I. 
Subsequently, we consider (5''uo)|K\(_;h,;fi) as originating from the finite volume 
discretization of two Cauchy-Dirichlet problems. For instance on (0,T) x (Z'\oo), 
both the initial condition Mo I (jh.oo) and the boundary condition Mh(-) := {S^uq){-,1'^) 
have the variation controlled in terms of Var uq and of j . It follows that the space- 
time variation of 5''wo on (0,r) x {R\{-l,l)) C (0,T) x (M \ Z'*)) is bounded 
uniformly in h; we refer to [ini Lemma 4.2], [51] Lemmas 5.3, 5.4] for the details. 
Using in addition the L°° bound, with the help of a diagonal extraction argument, 
we get convergence as ft. of a (not labelled) sequence {S'^uo)h to some limit u. 

Now we justify that the limit u is the unique fj-entropy solution of ( |1.7[ ), ( |1.2 ). 
The standard "weak BV" estimates for monotone finite volume schemes (see [33 ) 
permit us to get an approximate weak formulation and pass to the limit as /i — >■ 0. 



It follows that M is a weak solution to (|l.7|,(|l.2|). In the same way, the technique of 
3] allows us to get the Kru 
id {x > 0}. 

According to Definition 



[33| allows us to get the Kruzhkov entropy inequalities for u in the domains {x < 0} 
and {x > 0}. 

3.15| and Proposition |3.17[ what remains is to prove 



the entropy inequality (3.8) with an arbitrary non- negative test function ^ (not 
necessarily zero on the interface {x — 0}) and with any pair (c',c'') G G- Let us 



give a proof using (6.101. To this end, notice that a key feature of the scheme, 
thanks to the use the Godunov fiux of the exact fJ-Riemann solver at the interface, 
is that the discrete solutions it'' S'^uq take the values 

(y «")(•) = 5I„eN"o 

at the interface (strong one-sided traces at the interface {x = 0}), which satisfy 

((yu'')(t), (Yu'^m) e g, for a.e. t > 0. 
Because Q is an L^D germ, it follows that 

{{^^u'^){t),c') > q'- {{Yu'^){t),c^) , y{c\cn € G, (6.12) 
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for a.e. t > 0. Therefore with the same arguments as in [33], using (6.10) and using 



in addition the interface dissipation inequahty (6.121, we get the required entropy 
inequahties (3.8 1 with (c',c'') e Q and zero remainder term. 



Summarizing our findings, a ^-entropy solution has been constructed for any 
compactly supported initial function uq in BV{M.). It is easy to generalize this 
result so that it covers all uq G Indeed, because of t he Lip schitz continuity 

assumption on /'''", the contraction principle of Theorem 3.19| can be localized 
(as in the original work of Kruzhkov [50 ); using this contraction principle, by 
truncation and regularization of Uq we construct a strongly compact sequence of 
approximations . Then we pass to the limit in this sequence of approximations 



using the formulation of Definition 3.15 This concludes the proof of the theorem. 

□ 

6.4. On convergence of approximate solutions without BV estimates. In 

this section, we assume that the existence of a C/-entropy solution to the problem 



1/71), ( 1.2 ) is already known. For existence results, we refer, e.g., to Sections 6.1 6.2 



|6.3| and to many of the references at the end of the paper. Under this assumption, 
we have the following general convergence result. 

Theorem 6.5. Assume that the fluxes /'''" are merely continuous and the associated 
maximal L^D germ Q is chosen in such a way that there exists a solution to the 



problem (1.71,(1.2) for all bounded initial functions uq. 



Suppose that for any e > we are given a map : L°°(M) i— !■ L°°(M+ x M) with 
the following properties: 

(Bl) for each uq e L°°{R), the family {S'^uo)e>o is bounded in L°°(M+ x M); 

(B2) ifuo{x) = c'll{j.<o} +c''ll{j.>o} with {c\c'') £ Q, 

then S^Uq converges to Uq a.e. on M_|_ x M, as e — >■ 0,' 

(B3) for all uo,uo G L°°(R) and for all nonnegative ^ G P(M+ x M), there 
exists Ti = ri{uo,uo,^,e), with ri — > as e — > 0, such that the following 
approximate Kato inequality holds: 



\uq -uo|^(0, •) > -ri{uo,uo, (,,£). 



(B4) for all Uq G L°°{R) and for all ( G x M), there exists r2 = r2(uo,C,£) 

with r2 — > O as e — 0, such that the following approximate weak formulation 
holds: 

{S' un^t + fix, uo)^ + I uo^(0,-) -r-2(uo,^,£). 

Jr 

(B5) for all uq G L°°(M), for all nonnegative ^ G V{R+ x (M \ {a; = 0})), and 
for all fc G M, there exists r^ — r3(uo, fc, ^, e), with — > as e — > 0, such 
that the following approximate Kruzhkov formulation holds away from the 
interface {x — 0}; 

/ / {\S'uo - fc| 6 + q(x, S'uo, k)^ + f \uo- k\ e(0, •) > -rgK, fc, £)■ 

JR+ JR JR 

Then, as £ 0, S'^uq converges a.e. on x M to the unique Q entropy solution 



of (1.7), (1.2) 



The applications we have in mind are mainly related to numerical approxima- 
tions. Yet it should be noted that the standard vanishing viscosity approximations 



of Sections 6.1 fulfill assumptions (B1)-(B5). 

In this context, the property (B2) is somewhat delicate: it clearly holds for 
(c',c'') G QvY, then it can be extended to (c',c'') G Qvv, using the fact that Qvv 
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is the closure of Gyv (^^^ Proposition 5.1 1 and using the continuity with respect to 



the initial function uq of the vanishing viscosity solver . 

Thanks to Theorem |6.5[ from the existence result of Theorem |6.4| we deduce 

Corollary 6.6. // the vanishing viscosity germ Qvv is complete, then the following 



assumption can be dropped from Theorem 6.1 



/' and /' are not affine on any interval Id [0, 1]. 



(6.13) 



The germ Qyv is known to be complete in many situations (see [371 13H1 HH 1301 
ins [SI])- We expect that the vanishing viscosity germ Gyv is complete whenever 
are compactly supported, and also in the case there exist two sequences of 
elementary Gyy-entropy solutions (6fc(-))/c>i and (i?fc(-))fc>i such that as fc — >■ oo, 
bk — i> — oo and Bk oo. 

In general, assumption (B2) appears to be too restrictive. Certainly it holds for 
numerical schemes using either the exact Riemann solver on the interface {x = 0}, 
or the associated Godunov flux, but it may become difficult to justify in other 



situations. For instance, we cannot simply combine Theorems 6.4 and|6.5|in order 



to drop assumption (6.13) from the statement of Theorem 6.3 this is so because 
solely the connection solution A ll{2;<o} + B ll{a;>o} and the two constant states 
and 1 are preserved by the adapted viscosity approximation (6.6). However, the 



investigation of the approximate solutions for Riemann initial data of the form 
uo = c']l{,<o} +c' ll{^>oj, with ic\cn e ^(AB) \ {(AS), (0,0), (1,1)}, would 
require subtler arguments. 

Note that whenever the strong compactness property is available, the schemes 
that only preserve some definite germ Qq strictly smaller than Q do converge, thanks 



to the formulation of Proposition 3.17 (see [37], [S], [H] for a few particular cases). 



Proof of Theorem \6.5\ First, in view of (Bl), we have nonlinear weak-* compact- 
ness of the sequence {S^uq)^^^. Combining it with the strong convergence of S'^ug, 
with vP being any stationary solution from (B2), we pass to the limit as e — (up 
to extraction of a subsequence) in the formulations (B3), (B4), and (B5). Then the 
nonlinear weak-* limit fi of (a subsequence of) S^uq is a C/-entropy process solution, 
according to Remark |3.25| Applying Theorem |3.28 we conclude that /i is identified 
with the unique Q-entiopy solution u of (1.7), (1.2), and u is therefore the unique 
accumulation point of {S^Uo)^yo in Lj'Q^(M_|_ x M). This concludes the proof. □ 
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